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ABSTRACT 

The High-z Supernova Search Team has discovered and observed 8 new supernovae in the 
redshift interval z = 0.3-1.2. These independent observations, analyzed by similar but distinct 
methods, confirm the result of Riess et al. (1998a) and Perlmutter et al. (1999) that supernova 
luminosity distances imply an accelerating universe. More importantly, they extend the redshift 
range of consistently observed SN la to z rs 1, where the signature of cosmological effects has the 
opposite sign of some plausible systematic effects. Consequently, these measurements not only 
provide another quantitative confirmation of the importance of dark energy, but also constitute 
a powerful qualitative test for the cosmological origin of cosmic acceleration. We find a rate for 
SN la of (1.4 ± 0.5) x 10~ 4 h 3 Mpc~ 3 yr _1 at a mean redshift of 0.5. We present distances and 
host extinctions for 230 SN la. These place the following constraints on cosmological quantities: 
if the equation of state parameter of the dark energy is w — —1, then ffo^o = 0.96 ± 0.04, and 

— 1.40m = 0.35 ±0.14. Including the constraint of a flat Universe, we find Hm = 0.28 ±0.05, 
independent of any large-scale structure measurements. Adopting a prior based on the 2dF 
redshift survey constraint on Qm and assuming a flat universe, we find that the equation of state 
parameter of the dark energy lies in the range —1.48 < w < —0.72 at 95% confidence. If we 
further assume that w > —1, we obtain w < —0.73 at 95% confidence. These constraints are 
similar in precision and in value to recent results reported using the WMAP satellite, also in 
combination with the 2dF redshift survey. 

Subject headings: galaxies: distances and redshifts — cosmology: distance scale — supernovae: general 
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1. Introduction 

1.1. SN la and the Accelerating Universe 

Discovering Type la supernovae (SN la) with 
the intent of measuring the history of cosmic ex- 
pansion began in the 1980s with pioneering ef- 
forts by N0rgaard-Nielsen et al. (1989). Their 
techniques, extended by the Supernova Cosmol- 
ogy Project (Perlmutter et al. 1995) and by the 
High-z Supernova Search Team (HZT; Schmidt et 
al. 1998), started to produce interesting results 
once large-format CCDs were introduced on fast 
telescopes. Efforts to improve the use of SN la as 
standard candles by Phillips (1993), Hamuy et al. 
(1995), and Riess, Press & Kirshner (1996) meant 
that data from a modest number of these objects 
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at z i=a 0.5 should produce a significant measure- 
ment of cosmic deceleration. Early results by Perl- 
mutter et al. (1997) favored a large deceleration 
which they attributed to SIm near 1. But subse- 
quent analysis of an augmented sample (Perlmut- 
ter et al. 1998) and independent work by the HZT 
(Garnavich et al. 1998a) showed that the decel- 
eration was small, and far from consistent with 
Om = I- 

Both groups expanded their samples and both 
reached the surprising conclusion that cosmic ex- 
pansion is accelerating (Riess et al. 1998a; Perl- 
mutter et al. 1999). Cosmic acceleration requires 
the presence of a large, hitherto undetected com- 
ponent of the Universe with negative pressure: the 
signature of a cosmological constant or other form 
of "dark energy" If this inference is correct, it 
points to a major gap in current understanding 
of the fundamental physics of gravity (e.g., Car- 
roll 2001, Padmanabhan 2002). The consequences 
of these astronomical observations for theoretical 
physics are important and have led to a large body 
of work related to the cosmological constant and 
its variants. For astronomers, this wide interest 
creates the obligation to test and repeat each step 
of this investigation that concludes that an unex- 
plained energy is the principal component of the 
Universe. 

Formally, the statistical confidence in cosmic 
acceleration is high — the inference of dark energy 
is not likely to result simply from random sta- 
tistical fluctuations. But systematic effects that 
change with cosmic epoch could masquerade as 
acceleration (Drell, Lorcdo, & Wasscrman 2000; 
Rowan- Robinson 2002). In this paper we not only 
provide a statistically independent sample of well- 
measured supernovae, but through the design of 
the search and execution of the follow-up, we ex- 
pand the redshift range of supernova measure- 
ments to the region z « 1 . This increase in redshift 
range is important because plausible systematic 
effects that depend on cosmic epoch, such as the 
age of the stellar population, the ambient chem- 
ical abundances, and the path length through a 
hypothesized intergalactic absorption (Rana 1979, 
1980; Aguirre 1999a,b) would all increase with red- 
shift. But, for plausible values of Sl\ and Om, 0.7 
and 0.3 for example, the matter-dominated decel- 
eration era would lie just beyond z = 1. As a 
result, the sign of the observed effect on luminos- 
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ity distance would change: at z w 0.5 supernovae 
are dimmer relative to an empty universe because 
of recent cosmic acceleration, but at z w 1 the in- 
tegrated cosmological effect diminishes, while sys- 
tematic effects are expected to be larger (Schmidt 
et al. 1998). 

Confidence that the Universe is dominated by 
dark energy has been boosted by recent observa- 
tions of the power spectrum of fluctuations in the 
cosmic microwave background (de Bernardis et al. 
2002; Spergel et al. 2003). Since the CMB ob- 
servations strongly favor Cltotal = 1 to high pre- 
cision (±0.02), and direct measurements of CIm 
from galaxy clusters seem to lie around 0.3 (Pea- 
cock et al. 2001), mere subtraction shows there 
is a need for significant dark energy that is not 
clustered with galaxies. However, this argument 
should not be used as an excuse to avoid scrutiny 
of each step in the supernova analysis. Supernovae 
provide the only qualitative signature of the accel- 
eration itself, through the relation of luminosity 
distance with redshift, and most of that effect is 
produced in the recent past, from z = 1 to the 
present, and not at redshift 1100 where the im- 
print on the CMB is formed. This is why more 
supernova data, better supernova data, and super- 
nova data over a wider redshift range are needed 
to confirm cosmic acceleration. This paper is a 
step in that direction. 

The most likely contaminants of the cosmo- 
logical signal from SN la are luminosity evolu- 
tion, gray intergalactic dust, gravitational lens- 
ing, or selection biases (see Riess 2000, Filippcnko 
& Riess 2001, and Leibundgut 2001 for reviews). 
These have the potential to cause an apparent 
dimming of high-rcdshift SN la that could mimic 
the effects of dark energy. But each of these ef- 
fects would also leave clues that we can detect. 
By searching for and limiting the additional ob- 
servable effects wc can find out whether these po- 
tential problems are important. 

If luminosity evolution somehow made high- 
redshift SN la intrinsically dimmer than their local 
counterparts, supernova spectra, colors, rise times 
and light-curve shapes should show some concomi- 
tant effects that result from the different veloci- 
ties, temperatures, and abundances of the ejecta. 
Comparison of spectra between low-redshift and 
high-rcdshift SN la (Coil et al. 2000) yields no 
significant difference, but the precision is low and 



the predictions from theory (Hoflich, Wheeler, & 
Thielemann 1998) are not easily translated into a 
limit on possible variations in luminosity. Some 
troubling differences in the intrinsic colors of the 
high and low-redshift samples have been pointed 
out by Falco et al. (1999) and Leibundgut (2001). 
Larger, well-observed samples, including the one 
reported here, will show whether this effect is real. 

Gray dust that absorbs without producing as 
much reddening as Galactic dust could dim high- 
redshift SN la without leaving a measurable im- 
print on the observed colors. Riess et al. (1998a) 
argued that dust of this sort would need to dim 
distant supernovae by 0.25 mag at z w 0.5 in a 
matter-only universe. To produce the dimming 
we attribute to acceleration, dust would also in- 
crease the variance of the apparent SN la luminosi- 
ties more than is observed. However, such dust, 
if smoothly dispersed between the galaxies, could 
appear degenerate with cosmic acceleration (Rana 
1979, 1980; Aguirre 1999a,b). Aguirre's electro- 
magnetic calculations of the scattering properties 
for intergalactic dust show that it can produce 
less reddening than galactic dust, but it cannot 
be perfectly gray. Near-infrared observations of 
one SN la at « 0.5 (Riess et al. 2000) do not show 
the presence of this form of dust. Further obser- 
vations over a wide wavelength range have been 
obtained by our team to construct a more strin- 
gent limit; these will be reported in a future paper 
(Jha et al. 2003a). Recent work by Paerels et al. 
(2002), which failed to detect X-ray scattering off 
gray dust around a z = 4.3 quasar, seems to indi- 
cate that gray, smoothly distributed, intergalactic 
dust has a density which is too low by a factor of 
10 to account for the 0.25 mag dimming seen in 
the SN la Hubble diagram. 

Selection biases could alter the cosmological in- 
ferences derived from SN la if the properties of dis- 
tant supernovae are systematically different from 
those of the supernovae selected nearby. Simple 
luminosity bias is minor because the scatter in su- 
pernova luminosities, after correction for the light- 
curve shape, is so small (Schmidt et al. 1998; 
Perlmutter et al. 1999; Riess et al. 1998a). How- 
ever, the nearby sample of SN la currently spans a 
larger range of extinctions and intrinsic luminosi- 
ties than has been probed for SN la at z 0.5. 
Most searches to date have only selected the tip 
of the iceberg: most supernovae at high redshift 
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lie below the sensitivity limit. We assume that we 
are drawing from the same population as nearby, 
but it would be prudent to test this rigorously. 
To make the SN evidence for dark energy robust 
against sample selection effects, we need more sen- 
sitive searches for SN la at z w 0.5 that could de- 
tect intrinsically dim or extinguished SN la so we 
can verify that the distant supernovae have a sim- 
ilar range of extinctions and intrinsic luminosities 
as the nearby sample. This was one goal of the 
1999 search reported here. 

Extending the data set to higher rcdshift is a 
more ambitious and difficult way to test for the 
cosmological origin of the observed dimming ef- 
fect. Evolution or dust would most naturally lead 
to increased dimming at higher redshift, while cos- 
mic deceleration in the early matter-dominated 
era (z > 1.2) would imprint the opposite sign on 
luminosity distances. A search at higher rcdshift 
demand that we search to fainter flux limits in 
bands that are shifted to the red to detect rest- 
frame B and V. In this paper, we describe a sen- 
sitive search in the R and / bands carried out at 
the Canada-France-Hawaii Telescope (CFHT) and 
at the Cerro Tololo Inter- American Observatory 
(CTIO) in 1999. For the highest redshifts, we de- 
tect flux emitted in the ultraviolet at the source. 
Members of the HZT have also embarked on an ex- 
tensive study of the [/-band properties of nearby 
supernovae that will help with the interpretation 
of these high-redshift objects (Jha 2003c). 

Performing these tests requires searching for 
SN la with a deeper magnitude limit in redder 
bands than previous searches. This approach can 
sample the full range of extinctions and luminosi- 
ties atz« 0.5 and test for a turn-down in the Hub- 
ble diagram at z > 1. Thus far, direct measure- 
ment of deceleration at early epochs rests in obser- 
vations of one SN la, SN 1997ff, at the remarkably 
high redshift of <~ 1.7. These observations match 
best with a dark-energy source for the observed 
behavior of SN la (Riess et al. 2001). But this sin- 
gle object represents just one data point isolated 
from the body of SN la observations — and there 
is evidence that this object could be significantly 
magnified by gravitational lensing (Bemtez et al. 
2002). A continuous sample from z rts 0.5 through 
z w 1 out to z > 1.5 is required to make this cru- 
cial test convincing. The present paper represents 
a step in that direction by extending the sample 



to z w 1. Future discoveries of very high- z su- 
pernovae using the Advanced Camera for Surveys 
on the Hubble Space Telescope (HST), successfully 
installed in 2002, will help bridge the gap from the 
high-redshift end. 

Section 2 of this paper describes the search, 
shows spectra of the supernovae, and provides our 
photometric results. Section 3 gives an account 
of the analysis including K-corrections, fits to the 
light curves, and luminosity distances. In §4 we 
discuss the inferred cosmological parameters, and 
in §5 we discuss how these results can be used to 
test for systematic errors in assessing cosmic ac- 
celeration from SN la observations. Some novel 
features of the analysis are described more thor- 
oughly in the Appendix. 

2. Observations 
2.1. Search 

The SN la reported here were discovered at 
the CFHT using the CFH-12K camera and at the 
CTIO 4-m Blanco telescope using the CTIO Mo- 
saic camera. The CFH-12K camera provides 0'.'206 
pixels and a field of view of 0.33 deg 2 , and the 
CTIO mosaic has 0'.'270 pixels with a field of view 
of 0.38 deg 2 . 

We obtained templates in 1999 October and ob- 
tained subsequent images in November to find new 
objects, plausibly supernovae, with a rise time in 
the observer frame of ~ 1 month. On the nights of 
1999 October 3 and 1999 October 7 (UT dates are 
used throughout this paper), we obtained CFH- 
12K images of 15 fields (5 deg 2 ) in the I and R 
bands with median seeing of 0765 and 0'.'72, re- 
spectively. We integrated for a total of 60 min- 
utes on 6 of the /-band fields and 3 of the R- 
band fields, and 30 minutes on the rest. Each 
integration was of 10 minutes duration, with the 
frames dithered by small offsets to help in remov- 
ing cosmic rays and CCD defects. The photomet- 
ric zero-points were approximately / = 35.0 and 
R = 35.6 mag for the 60 minute exposures — in 
other words, a source of that magnitude would 
produce 1 e~ . This permits detection of stars with 
m = 24 mag at a signal-to-noise ratio (S/N) of 12 
and 16, respectively, or detection of a supernova 
with m = 24 mag at a S/N of 9 and 11 in a dif- 
ference search. A first-epoch search scheduled in 
October at CTIO was clouded out. 
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Despite a generous time allocation, we ran into 
difficulties at CFHT a month later because of 
weather. The first night (1999 November 2) had 
clouds and bad seeing (> 1"). We accumulated 
60 minutes in the I band for each of 8 fields, 
but only one of these fields contributed to the 
search, yielding one SN la: SN 1999ff. The second 
night (1999 November 3) had a mean extinction 
of about 0.5 mag from clouds, but the seeing was 
very good (median 0'.'59). We obtained reasonably 
good data in the / band for 8 fields, with integra- 
tions of 60 minutes for five of these, and 30 min- 
utes for three. Three nights later (1999 Novem- 
ber 6) we searched at CTIO, using i?-band ob- 
servations from CFHT as the first-epoch observa- 
tions. Although the CTIO data enabled us to con- 
firm the SN la found in the I band at CFHT, and 
produced several candidates that required spec- 
troscopic follow-up, we found no new supernovae 
from these observations. Spectra of several CTIO 
candidates revealed the disappointing fact that 
our attempt to compare November CTIO obser- 
vations with October CFHT templates produced 
subtly false candidates where a mismatch between 
the CFHT and CTIO filters caused the difference 
between the two to give the illusion of a new ob- 
jects appearing between October and November. 
In fact, these were M stars or emission-line galax- 
ies whose Ha emission line fell at the red end of 
the i?-band filter. Only the comparison of CFHT 
data with CFHT templates provided genuine de- 
tections of supernovae. 

Our reduction and search procedure, detailed 
by Schmidt et al. (1998), consists of bias subtrac- 
tion, flatficlding, subtraction of /-band fringes, 
masking bad columns, and finally combining 
dithered images to reject cosmic rays and mov- 
ing objects. In each field, we identified stars, 
performed an astrometric solution, and then 
remapped the image onto a common coordinate 
system. (This mapping is flux conserving but uses 
a Jacobian to restore the photometric accuracy 
that is destroyed by flatficlding.) For each pair 
of images of a field, from the first (October) and 
second (November) epoch, we then used the al- 
gorithms described by Alard & Lupton (1998) 
to convolve the image with the better seeing 
into agreement with the worse image. Next we 
matched the flux levels and subtracted. 

These difference images were then searched au- 



tomatically for residuals resembling the point- 
spread function (PSF), with possible SN la flagged 
for further inspection by humans. Simultaneously, 
we searched all image pairs by inspecting the dif- 
ference images by eye. We found little differ- 
ence in the detection efficiency between the semi- 
automatic approach and inspection by experienced 
observers. For all the SN la candidates, we in- 
spected each of the individual exposures to ensure 
that the putative SN la was not a cosmic ray, mov- 
ing object, or CCD flaw. All 37 final candidates 
were tabulated along with estimated / and R mag- 
nitudes (from the CTIO observations), including 
the properties of the possible host galaxy and dis- 
tance from the center of the host. This candidate 
list became our observing list for subsequent spec- 
troscopic investigation. 

Different exposure times in different fields, vari- 
able extinction due to clouds, changes in seeing, 
and sensitivity variations among the chips in the 
CFH-12K mosaic make it hard to quantify a sin- 
gle detection threshold for genuine events in this 
search. We believe that, overall, our search is 95% 
effective for objects brighter than / = 23.5 mag, 
and substantially better than 50% at I = 24 mag. 
In some favorable cases we found objects as faint 
as / = 24.5 mag. 

We began observing our 37 candidate objects at 
the Keck-II telescope on 1999 November 8. Our 
aim was to determine the nature of each object and 
to measure its redshift. Of these 37, two turned 
out to be active galactic nuclei at redshifts of 1.47 
and 1.67; one was an M star and two were galaxies 
with Ha emission admitted because of mismatched 
R filters; five had disappeared, suggesting either 
that the discovery was well past maximum or the 
detection was spurious; four were judged to be 
Type II supernovae from spectra or blue color; 
two were too bright to be of interest for our pur- 
poses, I :=y 21 mag in bright host galaxies, possibly 
nearby SN II, for which we did not spend time to 
get spectra; three were SN la with 0.3 < z < 0.7 
which we chose not to follow; seven were faint can- 
didates for which we did not have time to get spec- 
tra; and 11 were the candidate SN la which we 
chose to follow. These are shown in Figure 1. The 
selection of objects was based on color informa- 
tion (we chose objects with colors consistent with 
SN la with any amount of reddening, in the range 
0.1 < z < 1.5); the amount of variation of the 
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Fig. 1. — Host galaxies for the eleven supernovae. 
Each image is 20" on a side and taken from an av- 
erage of several I-band frames. A 2" radius circle 
marks the position of the supernova. North is at 
the top and east to the left in each image. 



SN, with preference to objects that were not seen 
in the previous epoch, but that had varied signif- 
icantly; the host-galaxy brightness at the super- 
nova position, avoiding supernova candidates on 
extremely bright backgrounds; and host bright- 
ness/size — avoiding supernovae in low redshift 
galaxies (z < 0.1). These selection effects are all 
undesirable, but with limited telescope time, they 
were necessary compromises to achieve the goals 
of the search. 

Table 1 lists the J2000 coordinates of the SN la, 
the Galactic extinction from Schlegel, Finkbeiner, 
& Davis (1998), the number of photometric obser- 
vations in various bandpasses accumulated from 
various telescopes, the modified Julian date of 
maximum light, and our nickname for the super- 
nova, assigned before there was a designation for 
each event in the IAU Circulars (see Tonry et al. 



1999). The last three objects on the list are real 
objects which we pursued, but based on our spec- 
tra and subsequent light curves, the evidence is too 
weak to consider them SN la. SN 1999fi (Boris) 
was very close to SN 1999fj, so it did not cost extra 
observing time to get photometry, but the other 
two (SN 1999fo and SN 1999fuj fooled us into 
making 22 photometric observations and expend- 
ing many hours of Keck time to get spectra that 
were, in the end, inconclusive. Although these two 
are real events, with a definite flux increase from 
October to November and a decrease in flux after 
November, we have not shown they are SN la, and 
they are not included in the analysis. 

2.2. Spectral Observations and Reduc- 
tions 

We obtained spectra of our SN la candidates 
with LRIS (Oke et al. 1995) on the Keck-II tele- 
scope during five nights between 1999 November 8 
and 14. The seeing varied from 0'.'8 to l'.'l during 
the first three nights and was ~0'/5 on the last two 
nights. We used a l'.'O slit for all of our observa- 
tions, except for SN 1999fv where we used a 077 
slit on one night. We used a 150 line mm^ 1 grat- 
ing for the first half of the run, then switched to 
a 400 line mm -1 grating for the second half in or- 
der to better remove night-sky lines. The resulting 
spectral resolution is ^20 A for the 150 line mm -1 
grating and ~8 A for the 400 line mm -1 grating. 
The pixel size for the 150 and 400 line mm -1 grat- 
ings is ^5 A pix -1 and ^2 A pix -1 , respectively. 
The total exposure times for each SN la are listed 
in Table 2. We moved the object along the slit be- 
tween integrations to reduce the effects of fringing. 
The slit was oriented either near the parallactic 
angle, or to include the nucleus of the host galaxy 
or a nearby bright star to provide quantifiable as- 
trometry along the slit and to define the trace of 
a source along the spectral direction. 

When using the 150 line mm^ 1 grating, inter- 
nal flatfield exposures and standard stars were ob- 
served both with and without an order-blocking 
filter to remove contamination from second-order 
light. For the high-z SN la observations no order- 
blocking filter was needed, as there is very lit- 
tle light from the object or from the sky be- 
low an observed wavelength of 4000 A. We used 
BD+17°4708 as a standard star for the first two 
nights and Fcige 34 (Massey et al. 1988) for 
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Table 1 
Fall 1999 Observations 



SN name 


KA (JzUUUj 


Dec (^Jzuuu ) 


hi\ri— V ) 


AT 

Nobs 


NUD max 


Nickname 


SN 1999fw 


23:31:53.03 


+00:09:32.3 


0.039 


25 


51482 


Nell 


SN 1999m 


02:27:58.33 


+00:39:36.8 


0.031 


15 


51488 


Fearless Leader 


SN 1999ff 


02:33:54.39 


+00:32:55.6 


0.026 


25 


51494 


Bashful 


SN 1999fn 


04:14:03.88 


+04:17:55.0 


0.326 


46 


51498 


Scooby Doo 


SN 1999fj 


02:28:23.72 


+00:39:09.6 


0.030 


20 


51475 


Natasha 


SN 1999fm 


02:30:35.62 


+01:09:43.3 


0.023 


21 


51488 


Fred 


SN 1999fk 


02:28:53.88 


+01:16:24.2 


0.030 


22 


51492 


Velma 


SN 1999fv 


23:30:35.80 


+00:16:40.0 


0.041 


8 


51457 


Dudley Doright 


SN 1999fo 


04:14:45.75 


+06:38:34.3 


0.291 


11 




Gertie 


SN 19996 


02:28:11.67 


+00:43:39.3 


0.031 


12 




Boris 


SN 1999fu 


23:29:48.21 


+00:08:27.0 


0.048 


11 




Alvin 



the last three nights. Standard CCD processing 
and optimal spectral extraction were done with 
IRAF 16 . We used our own IDL routines to cal- 
ibrate the wavelengths and fluxes of the spectra 
and to correct for telluric absorption bands. For 
SN la that were observed on more than one night, 
the data were binned to the larger pixel size if 
different gratings were used, scaled to a common 
flux level, and combined with weights depending 
on the exposure time to form a single spectrum. 

2.3. SN la Spectra and Redshifts 

The SN la spectra are shown in Figure 2 at 
their observed wavelengths. For each SN la we 
list in Table 2 the redshift, the supernova type, 
the number of nights observed with Keck, the to- 
tal exposure time, and an explanation of how the 
redshift was derived. For SN la with narrow emis- 
sion or absorption lines from the host-galaxy light, 
we determined the redshift by using the centroid of 
the line(s) with a Gaussian fit. For the three SN la 
spectra without narrow lines from the host galaxy, 
we use the broad SN la features to determine the 
redshift by cross-correlating the high- z spectrum 
with low-z SN la spectra obtained near maximum 
light. We also include galaxy and M-star spec- 

IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under coop- 
erative agreement with the National Science Foundation 



tra among our templates for cross-correlation with 
spectra without narrow emission lines. Errors on 
the redshifts are ± 0.001 (la) when based on a nar- 
row line and ±0.01 when based on broad SN la 
features. The redshift of SN 1999fv has a large 
uncertainty (1.17-1.22) due to the low S/N of the 
spectrum (see Coil et al. 2000 for further discus- 
sion). 

Figure 3 presents the high-z spectra smoothed 
with a Savitsky-Golay filter of width 100 A, sorted 
by redshift, with the lowest-z SN la at the top. 
This polynomial smoothing filter preserves line 
features better than boxcar smoothing, which 
damps out peaks and valleys in a spectrum. The 
width of the smoothing filter is apparent in the 
[O II] A3727 emission line seen in many of the 
spectra. We also plot two low-z SN la (SN 1989B, 
Wells et al. 1994; SN 1992A, Kirshner et al. 
1993) and a low-z SNIc [SN 19941 dereddened by 
E(B-V) = 0.45 mag; Filippenko et al. 1995] for 
comparison, all with ages before or at maximum 
light. Features which distinguish a SN la (Filip- 
penko 1997) are deep Ca II H&K absorption near 
3750 A, the Si II A4130 dip blueshifted to 4000 A, 
and Fe II A4555 and/or Mg II A4481, the com- 
bination of which create a distinct double-bump 
feature centered at 4000 A. For the lower- z spec- 
tra we have enough wavelength coverage to see 
6150 A, where the Si II A6355 absorption feature 
is prominent in SN la, but beyond z w 0.4 this 
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feature becomes difficult to detect. Early SNIc 
spectra look similar to SN la blueward of 5000 A, 
but lack the prominent Si II dip at 4000 A which 
leads to the double-bump feature seen only in 
SN la. 

SN 1999fw, 1999fh, 1999ff, 1999m, and 1999fj 
are all clearly SN la. SN 1999fm is a Type 1 SN, 
but we cannot distinguish from the spectrum alone 
whether it is a Type la or Ic. SN 1999fk shows a 
definite rise around 4000 A and has some hints of a 
double peak at this location, but is not definitively 
a SN la. SN 1999fv shows some bumps centered at 
4000 A, which, as discussed in Coil ct al. (2000), 
are plausibly SN la features at a redshift of ~ 1.2. 
All eight of these supernovae are used in this pa- 
per. 

SN 19996, 1999fu, and 1999fo do not show any 
SN features in their spectra. For each of these 
systems, we have a redshift; for SN 1999fi and 
1999fu this comes from a single line we assume to 
be [O II] A3727. For SN 1999fo the spectrum ap- 
pears to have Ca II H&K absorption at z = 1.07. 
As stated above, these three objects are not used 
in the analysis. 

2.4. Properties of the Host Galaxies 

Host galaxies and the environment around each 
supernova are shown in Figure 1. The images are 
the average of the /-band images with the best 
seeing. No host was detected for SN 1999fm. 
SN 1999fj appears in a small group of galaxies 
and its host is assumed to be the elongated galaxy 
just east of the supernova. Host magnitudes were 
measured from images taken after the supernova 
had faded. For the brightest hosts, the flux was 
summed out to an isophote of 26 mag/square- 
arcsec, but for the faint hosts the flux was summed 
in a 4" diameter aperture. The results, without 
correction for extinction, are given in Table 3 . 

Offsets between the supernovae and their hosts 
are also given in Table 3, as is the projected sepa- 
ration in kpc assuming a flat A-dominated cosmol- 
ogy with flM = 0.3. The distribution of projected 
separations for these new supernovae is narrower 
than that of earlier searches (Farrah et al. 2002), 
meaning our supernovae were discovered closer to 
the hosts. This may be due to the greater depth 
of this search and better image quality as it at- 
tempted to detect objects at z > 1. 
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Fig. 2. — High-z supernova spectra shown at their 
observed wavelengths, smoothed by a 32-pixel run- 
ning median and a Gaussian of 4-pixel sigma. 

2.5. Photometric Calibrations and Reduc- 
tions 

The photometric data gathered of the super- 
novae discovered at CFHT and CTIO came from 
a wide variety of instruments and telescopes in- 
cluding Keck LRIS, Keck ESI, Keck NIRC, VLT 
FORS, VLT ISAAC, UKIRT UFTI, UH 2.2-m, 
VATT, WIYN, and HST. We accumulated a to- 
tal of 216 observations of our 11 candidate SN la 
over the span of approximately one year, in the 
V, R, I, Z, and J bandpasses, totaling ~ 200 
hours of exposure time. The seeing of the com- 
bined, remapped images had a median of 0'.'83, 
with quartiles at 0772 and 1704. 

Calibration of these data turned out to be espe- 
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Table 2 

Fall 1999 Spectroscopic Data 



SN 


name 


z 


Type 


# nights 


Exp. (sec) 


z determination 


SN 


1999fw 


0.278 


la 


1 


1200 


broad SN features 


SN 


1999m 


0.369 


la 


1 


1800 


[0 II] A3727 in SN spectrum 


SN 


1999ff 


0.455 


la 


1 


2200 


Balmer lines in galaxy spectrum 


SN 


1999fn 


0.477 


la 


2 


5900 


[0 II] A3727 in SN spectrum 


SN 


1999fj 


0.816 


la 


1 


3600 


[0 II] A3727 in SN spectrum 


SN 


1999fm 


0.95 


la? 


2 


6600 


broad SN features 


SN 


1999fk 


1.057 


Ia? 


2 


7000 


[O II] A3727 in SN spectrum 


SN 


1999fv 


1.17-1.22 


Ia 


2 


10000 


broad SN features 


SN 


1999fo 


1.07 


? 


3 


14200 


Ca H&K dip in spectrum 


SN 


19996 


1.10 


? 


1 


5000 


[O II] A3727 in spectrum 


SN 


1999fu 


1.135 




2 


7200 


[O II] A3727 in spectrum 



Table 3 
Fall 1999 Host Galaxies 



SN name 


Host R 


Host I 


Host Z 


Offset (") 


Proj. Sep. (kpc) 


SN 1999fw 


21.1 


20.4 




0.51 


2.3 


SN 1999m 


21.2 


20.8 


20.5 


0.03 


0.2 


SN 1999ff 


19.9 


19.2 


19.1 


2.01 


12.6 


SN 1999fn 


24.5 


23.5 




0.31 


2.0 


SN 1999fj 


23.3 


22.3 


21.9 


0.89 


7.2 


SN 1999fm 




> 24.5 








SN 1999fk 




21.9 


23.7 


1.36 


11.9 


SN 1999fv 


> 26.0 


25.8 




0.30 


2.7 


SN 1999fo 




23.1 


23.2 


0.19 


1.7 


SN 19996 


24.8 


23.1 


22.9 


0.16 


1.4 


SN 1999fu 


24.2 


23.6 


23.5 


0.35 


3.1 



The offset is the angular separation between the SN and host in arcsec, and 
the projected separation is in kpc. SN 1999fm had no detectable host. Ellipsis 
indicate no image available in that band. 
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cially challenging, due to a variety of difficulties. 
The photometric calibration of stellar sequences 
near our SN la was intended to come from ob- 
servations with the UH 2.2-m and the CTIO 1.5- 
m telescopes. To this end, we observed Landolt 
(1992) standards in V, i?, /, and Z on a vari- 
ety of photometric nights, and also observed our 
SN la fields to establish the magnitudes of local 
photometric reference stars. We took a sequence 
of short exposures of a constant flatfield to estab- 
lish the shutter timing error, and checked the lin- 
earity of the CCD by exposures of different dura- 
tion. For standard-star reductions we summed the 
flux (with background subtraction) of each star 
through a 14" aperture, and calculated an atmo- 
spheric extinction term for each night. A typi- 
cal scatter in this fit to airmass and color was 
0.02 mag, and arises from the usual causes: sky er- 
rors, imperfect flatficlding, changes in atmospheric 
transparency, CCD non-linearity, shutter timing 
errors, and PSF or scattered-light variations. In 
the case of the UH 2.2-m data, there were some 
observations which deviated significantly from the 
root-mean-square (rms) error of the majority of 
data. This was eventually tracked down to an er- 
ror in dome control, causing observations to be vi- 
gnetted by the dome shutter. The UH 2.2-m data 
are therefore not completely reliable photometric 
calibrations. The CTIO 1.5-m data had no such 
problems, and were taken on absolutely photomet- 
ric nights. However, they do not cover all fields, 
and the overlap in brightness between data with 
high photometric accuracy from the 1.5-m and the 
unsaturated data from the larger telescopes was 
not large. 

To bridge these calibration difficulties, we com- 
bined two additional calibration sources, the Sloan 
Digital Sky Survey (SDSS) (Stoughton et al. 2002) 
and the 2001 HZT campaign (Barris et al. 2002). 
The 2001 HZT campaign took great pains to ob- 
tain very high accuracy photometry at the CTIO 
1.5-m telescope, which was reduced in a manner 
identical to the 1999 CTIO data, as described 
above. These observations of more than 200 stars 
were used to calibrate a catalog of R, I, and Z 
photometry of about 4000 stars near a = 02^ 28 m , 
S = +00°35', derived from CFHT+12K mosaic 
images. The R and I bandpasses of these data 
were transformed to the standard Kron-Cousins 
system (Kron & Smith 1951; Cousins 1976), us- 



ing the CTIO derived values, and the Z band was 
transformed to the CTIO system, as described be- 
low. 

This 2001 HZT catalog and all of our super- 
novae at RA = 2S h and Q2 h are encompassed by 
the SDSS catalog. Using the matched stars at 
02' 1 between our 2001 HZT photometric catalog 
and the SDSS (using PSF-derived magnitudes of 
the Sloan data — not the default available on the 
website), we derived the following transformations 
from the SDSS g', r', i', and z' magnitudes and 
our Vega-based magnitudes: 

V — g' = +0.172 - 0.745(g' - r'), 
R-r' = -0.122 - 0.316(r' - i'), 
I-i' = -0.423 - 0.378(i' - z'), 
Z-z' = -0.511 - 0.037(r' - i'), 
Z-z' = -0.508 - 0.267(r' - i') + 0.424(i' - z'). 

The two-color fit to Z — z 1 is slightly better, but 
obscures the fact that Z and z' are nearly identi- 
cal, apart from a zero-point offset. 

For all colors, the scatter in the transformation 
for individual stars (which are not dominated by 
shot noise like the faint SDSS stars, or by sat- 
uration like the CFHT stars of roughly the same 
magnitude) is approximately 0.02 mag rms. These 
transformations enable us to calibrate the R, 7, 
and Z 1999 CFHT survey fields via the SDSS 
early-release data. 

Since we did not observe in the V band during 
our Fall 2001 campaign, we needed extra steps to 
derive the calibration of the V observations for 
our 1999 data. We used Landolt (1992) stars in 
our fields and derived a satisfactory relation for 
< (R-I) < 1 mag: 

{V-R) = -0.024+ 1.104(i?-7). 

For redder stars, there are significant differences 
in the U-band colors for giants and dwarfs, which 
we sought to avoid. This relation allowed us to 
produce "V" magnitudes for the Fall 2001 catalog, 
using the R and I magnitudes. We then compared 
the stars that were also available from SDSS, and 
derived a g' to V transformation. While individual 
stars can have substantial errors via these trans- 
formations, in the mean of hundreds of stars that 
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we used these transformations are better than 0.01 
mag. 

The Z band was both important and problem- 
atic for us. At z > 1, the rest-frame B band starts 
to shift beyond the observer's I band, and the sys- 
tematics of rest-frame [/-band observations were 
not yet well established at the time of these ob- 
servations (but see Jha 2003c), so we needed a 
color for our SN la based on bandpasses redder 
than /. There is a variety of Z bands; ours is 
a Vega-normalized bandpass defined by the CTIO 
natural system (RG850, two aluminum reflections, 
atmosphere, and a SITc CCD red cutoff), and it 
is quite well described as a sum of two Gaussians, 
one centered at 8800 A with a full-width at half- 
maximum (FWHM) of 280 A and a height of 0.86, 
and the other at 9520 A with a FWHM of 480 A 
and a height of 0.41. The Z bandpass and this 
approximation are illustrated in Figure 4. Our Z- 
band observations were calibrated by observing a 
series of Landolt stars, whose magnitudes were de- 
rived by integrating their spectrophotometry (N. 
Suntzeff, in preparation) with this bandpass. This 
system is defined to have (V — Z) = mag for 
Vega. 

The SDSS observations covered all our fields ex- 
cept for those containing SN 1999fn and SN 1999fo. 
For all SN la, except these two, we extracted 
SDSS stars which fell within the supernova 
fields, and used the above relations to produce 
Johnson/Kron-Cousins magnitudes for each star. 
For SN 1999fn and SN 1999fo, we used our UH 
2.2-m photometric observations, verified with a 
few overlapping stars from the 1999 CTIO 1.5- 
m observations, to set the R and I magnitudes 
of a stellar sequence around the SN. The V, Z 
magnitudes were set using color relations derived 
above. As stated above, these transformations, 
while good in the mean, may not be appropriate 
for unusual stars (or quasars), and consequently, 
our derived magnitudes are based relative to the 
sequence of stars, rather than weighting any indi- 
vidual star too highly. The adopted standard-star 
sequences in the vicinity of each supernova are 
given in Table 4. This calibration procedure is far 
from ideal, but cross-checks suggest it is as good 
as the Sloan zero-point uncertainty, 0.04 mag. 

We calibrated the J-band data from VLT and 
UKIRT using observations of ten standard stars 
from Persson et al. (1998). The Keck data were 



calibrated by transferring the photometric zero- 
points from the UKIRT observations to a galaxy 
in the field of SN 1999m. 
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Table 4 
Local Standard Stars 
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+00:18 
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.01 





.07 


20 


.86 





.11 


19 


.09 





.05 


18.38 
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23: 


:30:46 


.6 


+00:20 


:25 


21 


.80 





.06 


20 


.60 





.09 


18 


.96 





.04 


18.31 





.06 


SN1999fh 


02: 


:28:04. 


.7 


+00:40 


:21 


19 


.32 





.02 


18 


.79 





.03 


18 


.27 





.03 


18.09 





.05 
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.03 
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.35 





.03 


18 


.79 
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02: 


:28:05. 


.7 


+00:41 


:12 


22 


.54 





.11 


21 


.38 





.17 


19 


.65 





.07 


18.96 





.10 




02: 


:28:06. 


.7 


+00:41 


:30 


21 


.87 





.07 


20 


.79 





.10 


19 


.48 





.07 


18.97 





.10 




02: 
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.6 


+00:35 


:46 


21 
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20 
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Fig. 4. — Z bandpass used for these observations. 
The dashed line is an approximation of the ac- 
tual transmission, made up of two Gaussians as 
described in the text. 



Fig. 3. — High-z SN spectra, shifted to rest wave- 
lengths, smoothed by a Savitsky-Golay filter with 
a width of 100 A and compared to two low-z 
SN la (SN 1992A and SN 1989B) and a low-z 
SNIc (SN 19941). The high-z SN spectra are shown 
sorted by redshift, with the lowest-z SN at the top. 
Based on these spectra, we cannot conclude that 
SN 1999fi, 1999fo, and 1999fu are SN la. 
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Table 4 — Continued 



SN name 



RA 



Dec 



V 



± 



R 



± 



± 



± 



SN1999fo 



SN1999fj 



SN1999fw 



SN1999fn 



SN1999fk 



02 


30 


41 


5 


+01 


08 


03 


21 


02 





03 


19 


90 


0.05 


18 


44 





03 


17 


86 


0.04 


02 


30 


2G 


3 


+01 


08 


59 


17 


56 





02 


17 


05 


0.02 


16 


57 





02 


16 


41 


0.02 


02 


30 


25 


6 


+01 


09 


21 


19 


62 





02 


18 


48 


0.02 


16 


80 





02 


16 


14 


0.02 


02 


30 


4G 


3 


+01 


09 


31 


20 


66 





03 


19 


58 


0.04 


18 


16 





03 


17 


61 


0.03 


02 


30 


24 


2 


+01 


10 


57 


20 


93 





03 


19 


89 


0.05 


18 


47 





03 


17 


91 


0.04 


02 


30 


37 


7 


+01 


12 


00 


19 


01 





02 


18 


03 


0.02 


16 


88 





02 


16 


44 


0.02 


04 


14 


4G 


4 


+06 


38 


46 


18 


87 





03 








17 


40 





03 


17 


28 


0.03 


04 


14 


40 


4 


+06 


39 


46 


19 


17 





03 








17 


74 





03 


17 


61 


0.03 


04 


14 


41 


4 


+06 


38 


40 


21 


09 





03 








18 


26 





03 


17 


91 


0.03 


04 


14 


49 


G 


+06 


37 


35 


19 


45 





02 








18 


30 





03 


18 


25 


0.03 


04 


14 


4G 


5 


+06 


38 


29 


20 


13 





02 








18 


81 





03 


18 


74 


0.03 


04 


14 


42 


8 


+06 


37 


52 


22 


04 





07 








19 


14 





03 


18 


77 


0.03 


04 


14 


43 


6 


+06 


37 


23 


21 


61 





05 








19 


08 





03 


18 


81 


0.03 


04 


14 


45 


G 


+06 


37 


47 


21 


60 





05 








19 


12 





03 


18 


86 


0.03 


02 


28 


13 


9 


+00 


39 


05 


19 


80 





02 


18 


77 


0.03 


17 


46 





02 


16 


94 


0.03 


02 


28 


28 


1 


+00 


38 


25 


20 


81 





03 


19 


71 


0.04 


18 


05 





03 


17 


40 


0.03 


02 


28 


17 


3 


+00 


39 


06 


21 


20 





04 


20 


02 


0.05 


18 


36 





03 


17 


70 


0.04 


02 


28 


2G 


5 


+00 


40 


53 


20 


86 





03 


19 


78 


0.04 


18 


33 





03 


17 


76 


0.04 


02 


28 


19 





+00 


3G 


59 


21 


83 





06 


20 


69 


0.09 


19 


05 





05 


18 


41 


0.06 


02 


28 


31 


7 


+00 


37 


01 


20 


79 





03 


19 


89 


0.05 


18 


89 





04 


18 


51 


0.07 


02 


28 


15 


7 


+00 


40 


27 


23 


91 





30 


22 


51 


0.48 


19 


86 





09 


18 


78 


0.08 


02 


28 


20 


3 


+00 


39 


40 


21 


72 





06 


20 


74 


0.10 


19 


46 





06 


18 


96 


0.10 


23 


31 


3G 


2 


+00 


15 


10 


21 


19 





04 


20 


00 


0.05 


18 


24 





03 


17 


54 


0.03 


23 


31 


37 


3 


+00 


10 


14 


18 


92 





02 


18 


50 


0.02 


18 


09 





03 


17 


95 


0.04 


23 


31 


5G 


4 


+00 


15 


47 


21 


19 





04 


20 


24 


0.06 


19 


18 





05 


18 


77 


0.08 


23 


32 


01 


7 


+00 


13 


41 


21 


61 





05 


20 


55 


0.08 


18 


95 





04 


18 


31 


0.06 


23 


31 


58 





+00 


14 


12 


20 


69 





03 


19 


67 


0.04 


18 


27 





03 


17 


72 


0.04 


23 


31 


58 


3 


+00 


13 


09 


20 


57 





03 


19 


51 


0.04 


18 


21 





03 


17 


71 


0.04 


23 


31 


52 


G 


+00 


09 


32 


20 


22 





03 


19 


56 


0.04 


18 


89 





04 


18 


66 


0.08 


23 


31 


59 


9 


+00 


11 


08 


20 


79 





03 


19 


66 


0.04 


18 


07 





03 


17 


45 


0.03 


04 


14 


09 


2 


+04 


19 


06 


18 


35 





03 


17 


79 


0.03 


17 


19 





03 


17 


19 


0.03 


04 


13 


59 


9 


+04 


18 


46 


19 


47 





03 


18 


24 


0.03 


16 


74 





03 


16 


44 


0.03 


04 


14 


07 


8 


+04 


19 


15 


19 


40 





03 


18 


24 


0.03 


17 


14 





03 


16 


91 


0.03 


04 


13 


59 


7 


+04 


19 


08 


19 


90 





05 


18 


68 


0.03 


17 


42 





03 


17 


11 


0.03 


04 


14 


07 


9 


+04 


17 


58 


19 


92 





02 


18 


75 


0.03 


17 


60 





03 


17 


37 


0.03 


04 


14 


05 


6 


+04 


17 


40 


20 


67 





03 


19 


41 


0.03 


18 


08 





03 


17 


79 


0.03 


04 


13 


58 


6 


+04 


1G 


37 


20 


81 





04 


19 


63 


0.03 


18 


20 





03 


17 


88 


0.03 


04 


14 


02 





+04 


17 


41 


21 


22 





05 


19 


94 


0.03 


18 


51 





03 


18 


20 


0.03 


02 


28 


48 


2 


+01 


1G 


17 


19 


33 





02 


18 


43 


0.02 


17 


43 





02 


17 


05 


0.03 


02 


28 


45 


G 


+01 


15 


39 


19 


96 





02 


18 


90 


0.03 


17 


47 





02 


16 


91 


0.03 


02 


29 


01 


1 


+01 


13 


42 


21 


89 





06 


20 


57 


0.08 


18 


59 





03 


17 


80 


0.04 


02 


28 


54 


G 


+01 


14 


43 


20 


79 





03 


19 


86 


0.05 


18 


81 





04 


18 


41 


0.06 
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All our supernova images were compared with 
the above stellar sequences to derive zero-points. 
The star fluxes were obtained using the Vista 
"psf" routine, which sums up the light within a 
radius of 20 pixels around a star, subtracts a fitted 
sky value, and rejects contamination from neigh- 
boring objects. Comparison of the fluxes and mag- 
nitudes, with due regard for saturation and mag- 
nitude errors, gave us a flux-magnitude zero-point 
and error for most images. For HST images, we 
used the precepts of Dolphin (2000) to convert 
fluxes to magnitudes. 

2.6. Photometry of SN la 

Our reductions of deep supernova observations 
obtained after the candidates were identified are 
very similar to the reductions carried out dur- 
ing the search. After bias subtraction, flatfield- 
ing, J-band and Z-band fringe correction, masking 
bad columns, and combination of dithered images 
while rejecting cosmic rays and moving objects, 
we identified stars relative to a fiducial image and 
performed an astrometric solution. We remapped 
a 1024 x 1024 pixel subarray of the source image 
onto a tangent plane coordinate system with 072 
pixels. The HST observations had a small enough 
field of view and large enough distortion that some 
amount of star selection by hand was necessary to 
obtain a satisfactory astrometric match. 

We had expected to use HST to acquire very 
high-quality light curves for a subset of our ob- 
jects, but the gyro failure in Fall 1999 made HST 
unavailable until after this group of SN la had dis- 
appeared. We did succeed in using HST to get 
template observations of the host galaxies after 
the supernovae had faded for most of our targets. 
Late-time templates are very helpful to subtract 
the host galaxy accurately and to set the SN la 
flux zero-point. 

Most modern supernova light curves are derived 
by subtraction of a template image taken before 
the supernova explosion or long after, so that it 
carries no supernova flux. While this sets a good 
flux zero-point for the subtraction, it also inflicts 
a common, systematic flux error on the light curve 
which can be substantial if the template observa- 
tion has poor S/N or very different seeing. 

To avoid these correlated errors, we subtracted 
all N(N — l)/2 pairs of observations from one an- 



other. We carried this out using the Alard & Lup- 
ton (1998) code mentioned previously. This pro- 
cedure produces negative fluxes for some subtrac- 
tions, but it has the advantage that every observa- 
tion serves as a template for the entire light curve. 
The details of this procedure are in the Appendix. 
We have found that this procedure reduces the er- 
rors by a factor of about \[2 (Novicki & Tonry 
2000). 

There are two unresolved questions after the 
N(N — l)/2 procedure has been followed. We do 
not know a flux zero-point, since it has been re- 
moved when we difference images, and we need 
a flux-to-magnitude conversion which incorpo- 
rates the independent magnitude zero-points es- 
tablished for each image. 

The flux zero-point is generally established by 
late-time observations when there is no remaining 
supernova flux. This is an improvement over the 
usual "template" method, because the flux zero- 
point which comes out of the N(N — l)/2 proce- 
dure for the late-time point comes from the av- 
erage comparison with all the other observations, 
and the relative fluxes of all the other observa- 
tions are not dependent solely on the final observa- 
tion. In some cases where we had several late-time 
points without supernova flux, we adjusted the 
flux zero-point to be an average of those points, 
or to include a first point that was obtained prior 
to the supernova explosion. 

Infrared template images were available only for 
SN 1999fm and SN 1999fk, so we used standard 
PSF photometry with the DaoPhot II/Allstar 
package (Stetson 1987; Stetson & Harris 1988) 
to measure the J-band magnitudes of the SN cor- 
rected to the same aperture used for the standard 
stars. SN 1999fm was not visible in the J band. 
SN 1999fk is near the edge of its host so contam- 
ination from the galaxy is negligible. SN 1999fn 
and SN 1999fv are embedded in their hosts, so 
contamination may be significant. SN 1999ff is lo- 
cated in a bright (J ss 17.5 mag) elliptical galaxy 
for which it was essential to subtract the host. 
We used the iraf / STSDAS analysis isophote tasks 
to construct a model of the host, subtracted our 
model, and then performed PSF photometry on 
the supernova. We were unable to define a PSF 
for the Keck observations of SN 1999fn due to the 
lack of suitable stars in the Keck/NIRC field of 
view, so we performed aperture photometry in an 
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Table 4 — Continued 



SN name RA 


Dec 


V 


± 


R 


± 


I 


± 


z 


± 


02:29:05.3 


+01:13:46 


21.50 


0.05 


20.46 


0.08 


18.95 


0.04 


18.35 


0.06 


02:28:42.9 


+01:15:31 


20.21 


0.03 


19.92 


0.05 


19.45 


0.06 


19.29 


0.13 


02:29:02.1 


+01:13:14 


22.63 


0.12 


21.49 


0.19 


19.92 


0.10 


19.30 


0.13 


02:28:42.3 


+01:15:16 


20.69 


0.04 


20.32 


0.07 


19.95 


0.10 


19.84 


0.22 



aperture with a radius of one arcsecond. 

Tables 5-12 give the derived supernova light 
curves, and Figure 5 shows the data points and 
typical fitted light curves. 
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Table 5 
Observations of SN 1999fw 



MJD 


m 




± 


K 


± 


obs 








V 




Kv^B 








51496. 


.25 


22. 


,18 


0.15 


-0.34 





.04 


991114_uh 


51514. 


.20 


23 


.84 


0.86 


-0.24 





.03 


991202_uh 


51529. 


.21 




. u 









.04 


QQ1 91 7 nVi 


51875. 


.25 


>24. 


.50 




-0.26 





.05 


001127_esi 








R 




Kr^v 








51455. 


.27 


>24. 


.50 




-0.20 





.15 


991004_12k 


51488. 


.15 


21. 


,30 


0.01 


-0.42 





.04 


991106_8k 


51496. 


.22 


21. 


,43 


0.02 


-0.30 





.06 


991114_uh 


51514. 


.26 


22. 


,40 


0.13 


-0.15 





.06 


991202_uh 


51517. 


.21 


22. 


66 


0.09 


-0.14 





.07 


991205Jris 


51529 


.28 


23 


,54 


0.21 


-0.09 





.07 


991217_uh 


51875. 


.26 


>24. 


,50 




-0.43 





.09 


001127_esi 








I 




Kj^R 








51454. 


.28 


>24. 


,50 




-0.53 





.05 


991003_12k 


51484. 


.34 


21. 


,16 


0.03 


-0.54 





.01 


991102_12k 


51485. 


33 


21. 


,14 


0.04 


-0.54 





.01 


991103_12k 


51496. 


.23 


21. 


,34 


0.12 


-0.52 





.03 


991114_uh 


51512. 


.31 


21. 


,98 


0.07 


-0.56 





.04 


991130_uh 


51517. 


.21 


22. 


.03 


0.03 


-0.57 





.04 


991205Jris 


51875. 


.28 


>24. 


,50 




-0.55 





.10 


001127_esi 








Z 




K z ^i 








51496. 


,28 


>24. 


,00 




-0.11 





.07 


991114_uh 


51517. 


.20 


21. 


,92 


0.48 


-0.08 





.15 


991205Jris 


51527 


.22 


>24. 


,50 




0.00 





.16 


991216_ufti 
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Table 6 
Observations of SN 1999fh 



MJD 


ra 




± 


K 




± 


obs 








R 




Kr^b 








51455. 


.47 


>24. 


.50 




-0. 


.85 





.12 


991004.12k 


51488, 


.25 


22 


,74 


0.02 


-0. 


.96 





.07 


991105.8k 


51496 


.51 


23 


.02 


0.10 


-1. 


.02 





.11 


991114.uh 


51528. 


23 


>24, 


,50 




-1 


,44 





.10 


991216.uh 


51570. 


.29 


>24. 


.50 




-1. 


.32 





.07 


000127Jris 


51875, 


.30 


>24, 


,50 




-1. 


,28 





.20 


001127.esi 








I 




Kl 


-^v 








51454. 


.51 


>24. 


,50 




-0, 


.88 





.07 


991003.12k 


51485, 


.36 


22. 


,40 


0.07 


-0. 


.92 





.02 


991103.12k 


51496 


.46 


22. 


,41 


0.07 


-0. 


,91 





.05 


991114.uh 


51570 


.27 


>24. 


.50 




-0. 


.94 





.07 


000127Jris 


51875, 


.29 


>24 


,50 




-0. 


.95 





.10 


001127.esi 
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Table 7 
Observations of SN 1999ff 



MJD m ± K ± obs 









R 






Kr^b 








51455. 


.41 


>24 


.50 






-0.68 


0. 


.05 


991004.12k 


51496 


.26 


22. 


.66 





.04 


-0.70 


0. 


.02 


991114.8k 


51496 


.35 


22 


,72 





.08 


-0.70 


0. 


.02 


991114_uh 


51497. 


.48 


22. 


.58 





.03 


-0.70 


0. 


.02 


991115_uh 


51517. 


.34 


24. 


.16 





.22 


-0.77 





.03 


991205Jris 


51529. 


.35 


24. 


.39 





.37 


-0.82 


0. 


.03 


991217_uh 


51530 


.28 


24. 


.83 





.52 


-0.82 


0. 


.03 


991218_uh 


51876. 


.26 


>24 


.50 
/ 






-0.78 
Ki^v 


0. 


.10 


001128_esi 


51454. 


.49 


>24 


.50 






-0.82 


0. 


.05 


991003.12k 


51484. 


.46 


22 


,86 





.07 


-0.84 


0. 


.00 


991102.12k 


51496. 


.28 


22 


.49 





.09 


-0.84 


0. 


.01 


991114.8k 


51496. 


.36 


22. 


.37 





.12 


-0.84 


0. 


.01 


991114.uh 


51497. 


.47 


22. 


,27 





.12 


-0.84 


0. 


.01 


991115.uh 


51517. 


.32 


23 


,26 





.09 


-0.84 


0. 


.01 


991205Jris 


51528. 


.30 


23. 


.59 





.20 


-0.85 


0. 


.01 


991216.uh 


51876. 


.27 


>24 


.50 

Z 






-0.86 

K Z ^R 


0. 


.05 


001128_esi 


51496 


.37 


22. 


.85 





.18 


-0.79 


0. 


.03 


991114.uh 


51517. 


.37 


22. 


.65 





.14 


-0.85 





.07 


991205_lris 


51527. 


.30 


23 


,51 


1 


.19 


-0.91 





06 


991216.ufti 


51873. 


.22 


23 


.67 
J 


1 


.56 


-0.94 
Kj^i 


0. 


.10 


001125.wiyn 


51501. 


.29 


22. 


.61 





.10 


-0.76 





.03 


991119_nirc 


51526. 


.31 


23 


,03 





.23 


-0.79 





.06 


991214_nirc 
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Table 8 
Observations of SN 1999fn 



MJD 


rn 




± 


K 




± 


obs 








73 
11 




K V, 


->B 








51459. 


.59 


^94 
^ZQ 


^.o 




n 

— U. 


7^ 
. 1 


0. 


.05 


QQ1 008 1 91r 
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3. SN la Analysis 
3.1. K-Corrections 

K-corrections were calculated using the for- 
mulas described by Kim, Goobar, & Perlmuttcr 
(1996) and Schmidt et al. (1998). The appar- 
ent brightness of a supernova observed in filter j, 
nij(t), but with a z — absolute magnitude light 
curve in filter i, Mi(t), is given by 



.//,•(/) = .U,-U1+ :)/')+2r,+.-,l„fv(^L) + A- ;y .((l + :)/'). 



where Kij is given by 



Kij = 2.5 log 



(1 + *) 



J F x (X)Si(X)dX 



/F A (A/(l + z))^(A)dA 



a) 



for filter energy sensitivity functions Si and Sj, 
with zero-points Zi and -2^, and supernova spec- 
trum F\. The luminosity distance, Dl, depends 
on the cosmological parameters, and is discussed 
extensively by Schmidt ct al. (1998). 

For multi-filter light curve shape (MLCS) fit- 
ting, K-corrections were calculated in the same 
manner as for Riess et al. (1998a), as part of 
the fitting process, using a prescription similar to 



that of Nugent et al. (2002). For K-corrections 
presented in Tables 5-12, we used a set of 135 
SN la spectra ranging from 14 d before maximum 
light to 92 d after maximum light. Extending 
the work of Nugent ct al. (2002), before apply- 
ing the K-correction formulae, the SN la spectra 
of each epoch were matched to the de-reddened 
U, B, V, R, I color curves of a well-observed set of 
SN templates, by warping the spectral energy re- 
sponse with a spline containing knots at the ef- 
fective wavelengths of each filter. In our case we 
used the compilation of Jha (2002) for SN 1991T, 
1994D, 1996X, 1997bp, 1997bq, 1997br, 1998ab, 
1998aq, 1998bu, 1998es, 1998V, 1999aa, 1999ac, 
1999by, 1999dq, 1999gp. The spectra were cor- 
rected for the estimated host extinction in the rest 
frame, as well as for Galactic extinction in the 
observer's frame using the Schlegel et al. (1998) 
values and the Savage & Mathis (1979) redden- 
ing law. The modified spectra were then used to 
calculate the K-corrections, providing a series of 
K-correction estimates as a function of SN age, 
and of SN template type. These values were fitted 
with a 3-knot spline, and uncertainties estimated 
from the scatter of allowable templates as dictated 
by the light-curve fits (in this case, the modified 
dml5 method of Germany et al. 2001). 

Our K-corrcctions were used to fit the light 
curve and update the time of maximum light, 
as well as the allowable supernova template type 
and host extinction, and the whole process was 
iterated until convergence (only two iterations 
were required). Tables 5-12 list the derived K- 
corrections for each observation along with their 
uncertainties. 



3.2. Light-Curve Fits using MLCS 

We fit the photometry of the eight SN la shown 
in Figure 5 using either the MLCS algorithm as 
described by Riess et al. (1998a) or Riess et al. 
(2001), depending on the quality of the data. 

For SN la whose light curve constrains the 
time of maximum well and whose colors provide 
sufficient leverage on potential reddening, we fit 
the photometry using the MLCS algorithm as de- 
scribed by Riess et al. (1998a). This method em- 
ploys the observed correlation between light-curve 
shape, luminosity, and colors of SN la in the Hub- 
ble flow, leading to significant improvements in the 
precision of distance estimates derived from SN la 
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Table 8 — Continued 
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± 
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Table 9 
Observations of SN 1999fj 



MJD m ± K ± obs 
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Table 10 
Observations of SN 1999fm 



MJD m ± K ± obs 
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Table 11 
Observations of SN 1999fk 
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Table 12 
Observations of SN 1999fv 
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0.05 


991121Jsaac 
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(Hamuy et al. 1995, 1996b; Riess, Press, & Kirsh- 
ner 1995, 1996a). 

The empirical model for a SN la light and color 
curve is described by four parameters: a date of 
maximum (t), a luminosity difference (at maxi- 
mum, A), an apparent distance modulus (/is), 
and an extinction (Ab)- Due to the redshifts of the 
supernova host galaxies we first correct the super- 
nova light curves for Galactic extinction (Schlegel 
et al. 1998), and then determined the host-galaxy 
extinction. 

Although the best extinction estimates and 
most precise distances arc found from SN la light 
curves which have extended temporal and spectral 
coverage, ideally to ^100 d after maximum and 
4 or 5 bandpasses (Jha 2002), we only use data 
within 40 d of maximum light and in rest-frame 
B and V to match the data which can be sampled 
at high redshift. This should only compromise our 
accuracy slightly, however. 

For SN 1999fv and SN 1999fh, failure of the 
HST gyros led to data with less than optimal cov- 
erage in time. For these two objects we used the 
snapshot distance method (Riess et al. 1998b) 
which employs the observed SN la spectrum to 
constrain the age. This approach is significantly 
less precise than the distance from a well-sampled 
light curve. 

3.3. Light Curve Fits using Ami§{B) and 
dml5 

We also use the Ami 5 method of Phillips et 
al. (1999) to estimate a distance to each of our 
SN la. Originally, the method characterized the 
light curve by a parameter tied to the £>-band 
photometry: the decline in magnitudes for a SN 
la in the first 15 days after B-band maximum. 
With the accumulation of a larger number of well- 
sampled light curves, it became possible to expand 
the method so that Amis(i?) is determined from 
the B, V, and /-band light curves. The Phillips 
et al. method then compares the light curve of a 
new object with a small number (six) of template 
light curves. 17 



17 The six templates are based on the light curves of seven 
objects: SN 1991T, 1991bg, 1992A, 1992al, 1992bc, and 
1992bo+1993H. The method now compares the light curve 
of a new object with 6 template light curves (Hamuy et al. 
1996b). 
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Fig. 5.— MLCS fits to our eight SN la. The open 
diamonds represent K-corrected ]/-band fits, while 
the fillcd-in circles represent K-corrected B-band 
fits, offset by +1 mag, except for SN 1999fh, which 
is V + 1 and R. SN 1999fo, 19996, and 1999fu 
have no light curves because we could not establish 
from the spectra that these were indeed SN la, so 
they were not followed. Redshifts are indicated in 
parentheses. 

In contrast to the MLCS method of Riess et al. 
(1996, 1998a), the Ami 5 method does not produce 
a continuum of "synthetic" templates. There is a 
small number of discrete templates. The redden- 
ing of a SN la is determined from the B — V and 
V — I colors at maximum, and from the B — V 
colors between 30 and 90 days after T^-band max- 
imum. It is assumed that unreddened SN la of all 
decline rates show uniform late-time color curves, 
as found by Lira (1995). Jha (2002, §3) found that 
the "Lira law" holds well for the 44 SN la stud- 
ied by him. Only exceptional cases like SN 2000cx 
seem to deviate greatly from the relation (Li et al. 
2001; Candia et al. 2003). 

Distances using the modified dml5 method are 
described in detail by Germany (2001). In sum- 
mary, well-observed SN la (approximately 15 ob- 
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jects) are fitted with a spline in BVRI and as- 
signed a value of Anii5(S) based on the method 
of Phillips (1993) and Phillips et al. (1999). These 
templates are extended to explosion times using 
the rise-time measurements of Riess et al. (1999b). 
SN la at cz > 4000 km s _1 with extinction values 
from Phillips et al. (1999) are fitted with these 
templates, with the best-fitting template used to 
assign these objects a dml5 value and maximum- 
light values B max , V max , R max , and I max . The 
Hubble velocity then provides an absolute magni- 
tude in each band, from which Mb,vm,i vs dml5 
relationships are derived. The values found from 
these relationships are applied to the set of well- 
observed templates to create a set of absolute mag- 
nitude light curves in BVRI. To measure the dis- 
tance to a new supernova, each template is fitted 
to the data, marginalizing over time of maximum 
and Ay, to find the distance modulus of the ob- 
ject. The distance is then estimated as the x 2 ~ 
weighted average of each template distance, with 
the distance uncertainty derived from the union of 
all templates and distances that have acceptable 
X 2 fits. 

3.4. Light Curve Fits using BATM 

The Bayesian Adapted Template Match (BATM) 
Method, described by Tonry (2003), uses a set of 
approximately 20 nearby supernova light curves 
which are well observed in multiple colors and 
which span the range of absolute luminosity to 
predict what they would look like at an arbitrary 
redshift with an arbitrary host extinction through 
the actual filters used for observation. Since these 
SN la do not in general have a significant amount 
of spectrophotometry, the BATM method uses 
approximately 100 spectral energy distributions 
(SEDs) of other SN la at a variety of SN la ages, 
and warps the SEDs into agreement with the pho- 
tometry. Given each of the template predictions 
and a light curve of a new SN la, the BATM 
method then calculates the likelihood distribution 
as a function of explosion time t , host extinction 
Ay, and distance d, and marginalizes over i to 
get a likelihood as a function of Ay and d. This is 
multiplied by a prior for Ay which is a Gaussian 
of a = 0.2 mag for negative Ay, a Gaussian of 
a = 1.0 for positive Ay, and a delta function at 
Ay = of equal weight. Unfortunately, for these 
observations (and probably most extant high red- 



shift SN la distances) there is strong covariance 
between d and Ay and the choice of prior makes 
a non-negligible difference in the final estimate for 
d. The BATM method then combines the prob- 
ability distributions for all of the template SN la 
and marginalizes over Ay to get a final estimate 
for the distance. 

The BATM method is still undergoing im- 
provement and is not currently "tuned" to min- 
imize scatter of the SN la in the Hubble flow. 
For example, one might expect that certain band- 
passes or time spans during a supernova's light 
curve might be better predictors of distance (the 
MLCS method is designed to optimize those 
choices). The BATM method simply fits all the 
data available in the observed bandpasses (avoid- 
ing the use of K-corrections) without regard for 
how a set of SN la in the nearby Hubble flow 
behave. As a result, the scatter of the BATM 
method for these Hubble-flow SN la is somewhat 
larger than that of MLCS (0.22 mag versus 0.18 
mag), but it has a different, and perhaps smaller, 
set of systematic errors. 

3.5. Fall 1999 Distances 

We list the distance estimates and uncertainties 
based on our light curves from MLCS, Anii 5 (B), 
dml5, and BATM in Table 13, along with the 
empty universe distance modulus corresponding to 
each redshift. The uncertainties are formal error 
estimates derived from each method, and we do 
not attempt here to evaluate whether these esti- 
mates are accurate. Below, in Table 15, we do 
pay attention to the self consistency of the various 
error estimates. This assumes a nominal value of 
i?o = 65 km s _1 Mpc -1 , though the cosmological 
inferences do not generally depend on the value of 
Hq. The estimates are generally consistent with 
each other, and we list a "final" distance which 
is the median of the estimates of each of the four 
methods. The uncertainty of this "final" distance 
is taken to be the median of the errors of the con- 
tributing methods, because we regard the process 
of taking a median of different analysis methods 
not so much sharpening our precision as rejecting 
non- Gaussian results. 
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Table 13 
Distance Moduli (h = 0.65) 



SN name 


MLCS 


Ami 5 


dml5 


BATM 


Final 


fl = 


SN 1999fw 


40.97 0.31 


40.89 0.15 


40.80 0.24 


40.99 0.42 


40.93 0.27 


40.82 


SN 1999m 






41.77 0.26 


41.46 0.47 


41.61 0.36 


41.52 


SN 1999ff 


42.38 0.29 


42.60 0.22 


42.35 0.43 


42.61 0.32 


42.49 0.30 


42.05 


SN 1999fn 


42.69 0.26 


42.36 0.16 


41.98 0.40 


42.19 0.23 


42.27 0.24 


42.18 


SN 1999Q 


44.03 0.22 


43.59 0.37 


43.76 0.35 


43.55 0.43 


43.67 0.36 


43.62 


SN 1999fm 


44.18 0.22 


44.06 0.23 


43.84 0.37 


43.77 0.35 


43.95 0.29 


44.05 


SN 1999fk 


<44.23 0.30 




44.26 0.39 


44.17 0.36 


44.23 0.37 


44.36 


SN 1999fv 


<44.59 0.39 






<44.19 0.45 <44.39 0.42 


44.74 



3.6. Rates 

From observations of just three SN la at z w 0.4 
Pain et al. (1996) calculated a rate of SN la 
which amounts to 34 SN la yr -1 deg -2 , with a 
la uncertainty of a factor of ~ 2, for objects 
found in the range of 21.3 < R < 22.3 mag. 
A more recent estimate by Pain et al. (2002) is 
(1.53 ± 0.3) x 10~ 4 h 3 Mpc~ 3 yr" 1 at a mean 
redshift of 0.55. Cappcllaro et al. (1999) re- 
port a nearby rate of SN la of 0.36 ± 0.11 h 2 
SNu (where 1 SNu = 10 10 L B& per century) or 
(0.79 ± 0.24) x 10~ 4 h 3 Mpc~ 3 yr^ 1 , using a local 
luminosity density of p = 2.2 x 10 8 h L B@ Mpc~ 3 . 
[There is significant uncertainty in the luminosity 
density however: Marzke et al. (1998) obtained 
1.7 x 10 8 , and Folkes et al. (1999) found 2.7 x 10 8 .] 
Pain et al. (2002) report a very modest increase in 
the rate of SN la with redshift, perhaps tracking 
the star formation rate which Wilson et al. (2002) 
estimate as being proportional to (1 + z) 17 in a 
flat universe with Qm = 0.3. 

Although our Fall 1999 observing was focused 
on finding candidates for detailed study, it still 
contains useful information on supernova rates. 
Two facts help establish rates from this search. 
First, our observations were very deep in the I 
band (since we were primarily searching for z > 1 
SN la) and second, because we also sought to find 
good z « 0.5 SN la for HST observations, we ob- 
tained spectra for every plausible candidate SN la 
with z < 0.7. While we did not obtain light curves 
for all these candidates, we do know which were 
Type la, and we know their rcdshifts. This pro- 



vides the raw material for a computation of the 
supernova rate. 

At the outset we need to make some reason- 
able assumption about the underlying luminosity 
distribution of SN la, calculate what fraction of 
these we would see in our search, and then de- 
rive the supernova rate. We assume a luminos- 
ity function for SN la at maximum which con- 
sists of three Gaussians: 64% "normal" SN la 
with mean absolute magnitude M B = —19.5 mag 
(h = i? /(100 km s" 1 Mpc" 1 ) = 0.65) and disper- 
sion a = 0.45 mag, 20% "SN 1991T" SN la with 
Mb = —19.8 mag and dispersion a = 0.30 mag, 
and 16% "SN 1991bg" SN la with M B = -18.0 
mag and dispersion a = 0.50 mag (Li et al. 
2001a). 

Table 14 gives values that a "normal" SN la 
might achieve at maximum, derived from the col- 
ors of SN 1995D at maximum and the spectral 
energy distribution of SN 1994S, with parenthe- 
ses indicating magnitudes that are less accurate 
because they involve an extrapolation beyond the 
data for either of these two objects. 

We work in a space of redshift z and m max , 
where m is the bandpass of interest (I band here). 
We use the magnitudes of SN 1995D at maximum 
and the SED of SN 1994S to convert a point in 
this space to a value for the S-band (rest frame) 
maximum and therefore a density from our lumi- 
nosity function. We multiply in a volume factor 
dV/dz/dSl and a factor of (l+z) -1 because we are 
interested in the observer's rate of observing SN la. 
All of this is done for a flat universe with (SIm, ^a) 
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Table 14 

Peak Magnitudes of Normal SN Ia 



z 


B 


V 


R 


I 


Z 


J 


0.05 


17.37 


17.41 


17.36 


17.82 


17.72 


(17.70) 


0.10 


18.94 


18.88 


18.80 


19.20 


19.24 


(19.30) 


0.20 


(20.70) 


20.29 


20.27 


20.36 


20.75 


(20.96) 


0.30 


(21.94) 


21.08 


21.07 


20.97 


21.28 


(21.63) 


0.40 


(23.07) 


21.89 


21.62 


21.55 


21.67 


21.97 


0.50 


(24.06) 


(22.59) 


22.10 


22.01 


21.96 


22.35 


0.60 


(24.90) 


(23.26) 


22.59 


22.41 


22.37 


22.79 


0.70 


(25.53) 


(23.97) 


23.06 


22.70 


22.70 


23.06 


0.80 


(26.00) 


(24.75) 


(23.54) 


22.91 


23.01 


23.11 


0.90 


(26.39) 


(25.46) 


(24.06) 


23.12 


23.17 


23.20 


1.00 


(26.75) 


(26.03) 


(24.61) 


23.48 


23.39 


23.25 


1.10 


(27.08) 


(26.45) 


(25.14) 


23.78 


23.50 


23.38 


1.20 


(27.39) 


(26.78) 


(25.66) 


(24.19) 


23.72 


23.59 


1.30 


(27.67) 


(27.07) 


(26.13) 


(24.59) 


24.11 


23.72 


1.40 


(27.95) 


(27.35) 


(26.55) 


(24.91) 


(24.39) 


23.94 


1.50 


(28.20) 


(27.60) 


(26.93) 


(25.36) 


(24.70) 


24.08 


1.60 


(28.45) 


(27.85) 


(27.24) 


(25.85) 


(25.04) 


24.18 


1.70 


(28.68) 


(28.08) 


(27.51) 


(26.30) 


(25.32) 


24.32 


1.80 


(28.90) 


(28.30) 


(27.75) 


(26.77) 


(25.72) 


24.38 
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= (0.3 , 0.7). We then convolve the resulting 
distribution with a dust extinction distribution 
which approximates the detailed curves of Hatano, 
Branch, & Deaton (1998): 25% of hosts are "bulge 
systems" with the probability of a given host 
A B given by f(A B ) oc 0.02 8{A B ) + lO" 125 "^, 
75% of hosts are "disk systems" with probability 
f(A B ) oc 0.02 S(A B ) + e - 3 - 77 - A B . (The redshifted 
SEDs tell us how the observed /-band extinction 
depends on a given host galaxy extinction A B .) 
This gives us a distribution for the observed rate 
of SN la as a function of peak apparent /-band 
magnitude. 

However, two complications arise. The "obser- 
vation time" which must be multiplied into this 
distribution is a convolution of the search time and 
the intrinsic time that a given SN la is above the 
detection threshold. Also, our supernova detection 
threshold is not based only on peak /-band mag- 
nitude. We have two search epochs, separated by 
31 days, and we have a sensitivity of / ss 24 mag 
to the flux difference between the first and second 
epochs. 

Using the formalism of the BATM method, 
we place SN 1995D (Am 15 = 0.99, MLCS A = 
-0.42) and SN 1999by {Am 15 = 1.93, MLCS 
A = +1.48) in the z-m max space. Any point in 
the space between them is assigned a light curve 
which is the flux interpolation between the red- 
shifted light curves of SN 1995D and SN 1999by. 
Points which lie brighter or fainter than these two 
SN la are assigned copies of one or the other light 
curves which are simply scaled brighter or fainter 
to match the desired peak brightness. Given a 
light curve, we calculate the visibility time for each 
point in the (z, mi) diagram, using a 31 d lag time 
and a sensitivity of / < 24 mag. (Our actual prob- 
ability of detecting a supernova rolls off gradually 
at fainter magnitudes, but we had S/N w 10 at 
/ < 24 mag and we ignored candidates which we 
thought would peak at / > 24 mag, hence the 
approximation of a step at / = 24 mag.) 

The resulting difference light curve is identical 
to the actual light curve when the explosion oc- 
curs after the first epoch, but for explosions that 
erupted earlier than the first epoch, taking the dif- 
ference between our two observation epochs de- 
creases the difference and can be negative. Low- 
luminosity SN la with rapid rise and declines 
would be visible for a shorter time than higher- 



luminosity SN la with their slowly declining light 
curves, so they are less likely to be found. For the 
purpose of calculating rates, the total time that 
a supernova could be detected is a convolution of 
the time that the difference flux from the super- 
nova explosion exceeds the sensitivity limit and 
the search time (a delta function in this case) . 

Figure 6 shows contours of predicted rates 
(number per observer's month per square degree 
per magnitude per unit rcdshift). We plot our 
eight SN la as well as three other SN la which 
we did not follow. Without light curves for these 
three, we can only show the discovery (difference) 
magnitude which is an upper limit to the peak / 
band magnitude that was eventually reached. 

The normalization of Figure 6 is adjusted to 
bring the cumulative distribution as a function of 
z into agreement with our observations. Figure 7 
integrates the lower panel of Figure 6 over magni- 
tude and then assembles the cumulative predicted 
number of SN la as a function of redshift. We 
compare this with the number observed to normal- 
ize the model. Our sensitivity limit was approxi- 
mately / < 24 mag, and the time between our first 
and second search epochs was 31 d. Our total area 
surveyed (including the cloudy second epoch) was 
3.0 deg 2 , of which we lost about 15% due to chip 
edges, bad columns, stars, etc., so our effective 
search area was about 2.5 deg 2 . Subsequent tests 
of our search technique using synthetic SN la in- 
jected into the images indicate that we miss very 
few objects down to the quoted sensitivity limit. 

Uncertainties in the searched area, effective 
time sampling, and modeling of supernova prop- 
erties amount to perhaps 10-20%, and a similar 
uncertainty is the fraction of objects which we de- 
tected photometrically but for which we did not 
obtain spectroscopy. Since our goal was to find 
suitable targets for HST studies and to find su- 
pernovae at z w 0.5, we were not complete in our 
spectroscopic follow-up. Subjective factors go into 
the decision of which objects to follow. For exam- 
ple, if we find a candidate at / = 23.5 mag, but 
it appears to be in a galaxy at z w 0.3, we might 
very well choose not to get a spectrum, reason- 
ing that "the supernova is well past maximum" or 
"it might be a type II" or "it might be highly ex- 
tinguished, or a SN 1991bg-type" . Spectroscopic 
resources are so valuable that our first priority, to 
obtain targets, often outweighs our desire to con- 
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trol systcmatics. The factors in these decisions 
change with redshift of the supernova: we believe 
we were complete in our search to z < 0.6 (even 
though we did not follow three SN la photomet- 
rically), but we were definitely not complete for 
more distant and fainter SN la. The lower panel 
of Figure 6 shows that the counts of SN la to a 
redshift of z < 0.6 are not very sensitive to our 
estimate of our completeness limit. 

Figure 7 shows the comparison between this 
model cumulative distribution and the observed 
rates for our Fall 1999 campaign. The Poisson 
statistics of the 8 SN la counted with z < 0.6 
are still the largest source of uncertainty, and we 
estimate our overall normalization is accurate to 
about 35%. The slope agreement between the ob- 
servations and the model in Figure 7 reassures us 
that our model is reasonably accurate. 

These model predictions, which we calculate for 
comparison with the observations, include com- 
plicated factors for visibility within a time win- 
dow and sensitivity limit, convolution with dust 
extinction, time dilation factors of (1 + z), and 
volume factors of dV/dz/dQ. The fundamental 
normalization that goes into the calculation and 
that brings the model into agreement with the ob- 
servations shown in Figure 7 amounts to 40,000 
SN la of all sub-types per Gpc 3 per (rest frame) 
year for the Hubble constant used above (Ho = 
65 km s _1 Mpc^ 1 ) for the SN la luminosity func- 
tion. Put in more sensible units, this is a rate of 
r = (1.4 ± 0.5) x 10" 4 h 3 Mpc~ 3 yr" 1 SN la at a 
mean redshift of z = 0.46. So the observed rate of 
SN la from a shell at redshift z is 

dN = rdn (dV/dz/dfl) (1 + z)- 1 At, 

where At has to be integrated over the luminos- 
ity function of SN la and the search parameters. 
These rates are in excellent agreement with those 
of Pain et al. (2002), and are not inconsistent with 
the local rates from Cappellaro et al. (1999), par- 
ticularly given the uncertainty in the local lumi- 
nosity density. 

Although we were not complete in our counts 
of SN la at z ~ 1, we do not believe it is possible 
that the rate of SN la closely tracks the star for- 
mation rate (SFR). Application of the SFR from 
Wilson et al. (2002) would suggest that the rate 
at z w 1.1 should be three times as great as oc- 
curs locally and nearly twice as great as the rate 



at z = 0.46. Our model illustrated in Figure 7, 
if altered from constant rate per comoving vol- 
ume to one which tracks the SFR of Wilson et 
al. (2002), would predict 16 SN la deg~ 2 in our 
search instead of 6.5 (we discovered 4). It is our 
impression that the constant rate per volume is 
closer to the truth — SN la at mj = 24 mag are 
simply not that common. Obviously, SN la explo- 
sions lag star formation by some unknown time 
distribution. This suggests the lag must span a 
substantial fraction of the time between the peak 
of star formation and the time when a supernova 
observed at z ~ 1 explodes, which is about 2 Gyr. 
As pointed out by Pain et al. (2002), we should 
obtain better constraints on any increase in SN la 
rate at z > 1 from ongoing search campaigns. 

3.7. A Homogeneous Set of SN la Dis- 
tances 

As standardizable candles, SN la have played 
an extremely important role in recent years. How- 
ever, observations of SN la are so difficult and 
resource-intensive that the results have appeared 
in a large number of papers, using a variety of dis- 
tance zero-points, photometric calibrations, etc. 
We have therefore undertaken to draw all the re- 
sults available to us together into a single table 
with a common calibration. 

The most recent data apart from the obser- 
vations presented here are the results from the 
campaigns of the HZT in 1997 (Leibundgut et al. 
2003), 1998 (Suntzeff et al. 2003), 1999 (Clocchi- 
atti et al. 2003), and 2000 (Jha et al. 2003b). 

In the 1997 data, there are four spectroscopi- 
cally classified SN la (SN 1997as, 1997bb, 1997bh, 
and 1997bj) and one object in an elliptical galaxy 
(SN 1997bd) at redshifts in the range 0.33 to 0.67. 
The redshifts are all from narrow emission lines as- 
sociated with the host galaxies, and on the promi- 
nent Ca II H&K absorption in the elliptical galaxy. 
The supernova spectrum in all cases yielded an age 
concordant with the light curves. Light curves in 
the R and / filters have been established for all ob- 
jects. The campaign in 1997 made use of ground- 
based telescopes only. The distance moduli, col- 
ors, and absorptions are given in Table 15. 

From our 1998 campaign, Suntzeff et al. (2003) 
report photometry and spectroscopy of six SN la: 
19981, 1998J, 1998M, 1998ac, 1998ai, and 1998aj. 
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A seventh object (SN 1998ah) was only observed 
spectroscopically. These seven objects range in 
redshift from 0.43 to 0.887. One of the significant 
results of the 1998 and 1999 (see below) campaigns 
was the comparison of the photometry of the SN la 
and nearby field stars from ground-based and HST 
observations: we found that the zero-points of the 
photometry agree to within 0.02 mag, which is 
reassuring. Since the highest-weighted points in 
light curves of high-z supernovae are usually the 
HST data, owing to their smaller statistical un- 
certainties, we naturally wanted to compare pho- 
tometry obtained on the ground and in space. Be- 
cause the evidence from supernova observations 
for a non-zero cosmological constant hinges on a 
systematic faintness of SN la at z w 0.5 by 0.22 
mag, we are now quite confident that the system- 
atic faintness of supernovae at such redshifts is not 
due to a zero-point error between space-based and 
ground-based light curves. 

In the Spring 1999 campaign, Clocchiatti et al. 
(2003) find five spectroscopically classified SN la: 
SN 1999M, 1999N, 1999Q, 1999S, and 1999U. All 
of them were followed up from ground-based ob- 
servatories and two of of them, SN 1999Q and 
SN 1999U, were observed with HST. We report 
in Table 15 the distance moduli to the latter two 
based on the HST photometry. The redshift of 
SN 1999U was obtained from emission lines in the 
parent galaxy, and that of SN 1999Q from cross 
correlating its spectrum with nearby SN la close 
to maximum light. 

The 2000 campaign was intended primarily to 
investigate the possibility that gray dust was caus- 
ing the dimming we attribute to cosmology. Jha et 
al. (2003b) present photometry over a very wide 
wavelength range. There are six SN la at a redshift 
of ~ 0.5 whose distances we include in Table 15: 
SN 2000dz, 2000ea, 2000ec, 2000ee, 2000eg, and 
2000eh. 

In addition to these observing campaigns of the 
HZT, we also drew from the literature all SN la 
results we could find. In cases where we had access 
to the observational data and photometry we re- 
computed fits to the distance and extinction using 
the MLCS98 method of Riess et al. (1998a), the 
MLCS2k2 improvement of Jha et al. (2003c), and 
the dml5 and BATM methods described above. 
In other cases, such as the Supernova Cosmology 
Project (SCP; Perlmutter et al. 1999), which has 



not published their photometry, we took their pub- 
lished distances, extinctions, and uncertainties at 
face value. One might expect that applying mul- 
tiple fitting techniques as we have done helps de- 
crease errors in the derived distances. The SCP- 
only results appear to have at least a factor of 
2 more scatter with respect to a cosmological fit 
than our data. This could be due to the quality 
of the underlying photometry, or to a less varied 
approach to the analysis. 

Assembling all of these results from the liter- 
ature and performing new fits when possible, we 
obtain long lists of SN la results for each method 
which arc internally self-consistent, but which 
have unknown offsets with respect to one another 
because of differing and inconsistent assumptions 
about Hubble constants and zero-points. We also 
have one more distance measure, derived from the 
redshift of the SN la in the "Hubble flow" (defined 
as 0.01 < z < 0.1) and the luminosity distance for 
an empty universe: 

All of the methods overlap for many SN la, par- 
ticularly the Calan-Tololo sample (Hamuy et al. 
1996a), with even the SCP results having 17 or 
18 overlaps with the rest, and the other methods 
having 52-108 overlaps. In order to determine the 
zero-point offsets we compute the median differ- 
ence between all pairs of methods, forming an an- 
tisymmetric matrix, and then fit this matrix as dif- 
ferences between components of a zero-point vec- 
tor. The residuals are quite well behaved, with 
the worst agreement between SCP and MLCS98 
at 0.07 mag, and the rest 0.02 mag or better. We 
put all the methods on the Hubble flow zero-point, 
i.e., distance in Mpc multiplied by Ho- 
lt is beyond the scope of this paper to ana- 
lyze and optimize the combination of the pub- 
lished data. We do note that there actually does 
appear to be a modest improvement in accuracy 
(judged by scatter with respect to the cosmolog- 
ical line) by combining different fitting methods. 
Nevertheless, we are concerned about the suscep- 
tibility to bias which might creep in if we optimize 
scatter with respect to the parameters we seek to 
measure, so we choose to use medians, a simple 
and robust combination scheme. Our "best dis- 
tance" for each SN la is just the median of all the 
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distances of all the different methods (excepting 
redshift-based luminosity distance!). The best dis- 
tances for the Fall 1999 data presented here differ 
slightly from those in Table 13 because we did not 
have Amis distances for more than the Fall 1999 
data and we thought that for the sake of consis- 
tency we should treat the Fall 1999 data in the 
same way as the other SN la for which we have 
observations. The scatter between the distances 
in Table 13 and Table 15 is 0.08 mag. 

Obtaining accurate uncertainties is somewhat 
more complicated. The scatter of any method's re- 
sults around any cosmological may not match the 
error estimates. We therefore examine this scat- 
ter for each method, prune at 3cr, and multiply 
the error estimates for each method by a factor to 
bring them into agreement with the scatter. This 
appears to be a more realistic adjustment than 
adding a "cosmic scatter" in quadrature (which 
we already do anyway when called for). We then 
obtain a final error estimate as the median of all 
of the contributing methods' errors diminished by 
the 1/4 power of the number of contributors (a 
factor empirically derived by examining scatter 
around cosmological models). As above, we be- 
lieve the covariance between the different methods 
applied to the same data is such that we are get- 
ting only a minimal improvement in accuracy by 
taking a median, but we hope that we are thereby 
"Gaussianizing" the distribution. The final list 
comprises 230 SN la. 

We will use this list below for determining cos- 
mological parameters, and although we will not 
repeat the data we found in the literature, we 
thought it important to publish the actual dis- 
tances which are being used for calculation of cos- 
mological parameters as well as other relevant in- 
formation which may be of use to others. In Ta- 
ble 15 we list the supernova name, its Galactic co- 
ordinates, host, and literature reference. We give 
the redshift as log(cz) and the luminosity distance 
and error in the same units, i.e., we do not want 
to use uncertain units of Mpc when we know the 
Hubble flow zero-point so exquisitely well. The 
median of the estimates for host-galaxy extinc- 
tion in the rest-frame V band follows, where we 
have used a standard Cardelli, Clayton, & Mathis 
(1989) extinction law to bring estimates in differ- 
ent bandpasses to V. Finally we identify the meth- 
ods which contributed to the median: MLCS98 = 



R, SCP/stretch = P, MLCS2k2 = J, dml5 = S, 
and BATM = T. 
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Table 15 
SNIa Summary 



SN 


l" 


b n 




z 


Host 


Ref. a 


log( 


cz) 


(log(dflo)} 




± 


(Av) 


Contrib 


sn72E 


314 


.840 


30 


.080 





.0023 


N5253 


16 


2 


.839 


2.399 





.033 


0. 


,10 


RT 


sn80N 


240 


.161 


-56 


.689 





.0056 


N1316 


9 


3 


.225 


3.140 





.043 


0. 


.21 


RT 


sn81B 


292 


.970 


64 


.743 





.0072 


N4536 


2 


3 


.334 


3.077 





.041 





.25 


RT 


sn81D 


240 


.161 


-56 


.689 





.0056 


N1316 


9 


3 


.225 


3.044 





.055 


0. 


.31 


JT 


sn86G 


309 


.543 


19 


.401 





.0027 


N5128 


26 


2 


.908 


2.440 





.035 


0. 


,49 


RT 


sn88U 


8. 


.737 


-81 


.227 





.3100 


Anon 


24 


4 


.968 


5.096 





.072 


0. 


10 


RT 


sn89B 


241 


.991 


64 


.403 





.0036 


N3627 


37 


3 


.033 


2.844 





.030 


0. 


,74 


RST 


sn90N 


294 


.369 


75 


.987 





.0044 


N4639 


21 


3 


.120 


3.204 





.035 


0. 


,30 


RST 


sn90O 


37. 


.654 


28 


.360 





.0307 


M+034403 


10 


3 


.964 


3.977 





.025 


0. 


,09 


RPJST 


sn90T 


341 


.503 


-31 


.526 





.0400 


P63925 


10 


4 


.079 


4.101 





.042 


0. 


,25 


RJT 


sn90Y 


232 


.645 


-53 


.854 





.0390 


Anon 


10 


4 


.068 


3.985 





.039 


0. 


,54 


RJT 


sn90af 


330 


.823 


-42 


.235 





.0500 


Anon 


10 


4 


.176 


4.149 





.026 


0. 


.06 


RPJST 


sn91M 


30 


.392 


45 


.900 





.0076 


11151 


5 


3 


.358 


3.397 





.030 


0. 


.00 


R 


sn91S 


214 


.061 


57 


.426 





.0560 


U05691 


10 


4 


.225 


4.283 





.042 


0. 


.11 


RJT 


sn91T 


292 


.609 


65 


.191 





.0070 


N4527 


21 


3 


.322 


2.961 





.028 


0. 


.50 


RST 


sn91U 


311 


.824 


36 


,212 





.0331 


14232 


10 


3 


.997 
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.023 


0. 


.08 


RST 


sn92G 


184 


.623 


59 
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.030 


0. 
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RT 
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23 
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10 
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0. 


,27 


RJST 
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16 
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10 


3 
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.036 


0. 
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73 
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10 


3 
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0. 
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.706 


-41 
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10 


4 
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0. 
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38 
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10 
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0, 
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.346 
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10 
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.030 


0. 
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.0450 
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10 
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.026 


0. 
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.4580 
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5 
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.0580 
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10 


4 
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0. 


,11 
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10 


4 
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10 
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.021 


0. 
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RPJST 
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10 
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.026 


0. 
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10 


4 


.421 


4.441 
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0. 
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.028 


-55 


.345 





.0630 
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10 


4 
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4.327 





.026 


0. 


.19 
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273 


.324 


20 


.460 
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10 


4 
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4.345 





.026 


0, 


.32 
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sn93H 


318 


.223 


30 
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E445-66 


10 


3 
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.026 


0. 


.32 


RJST 
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312 


.414 


28 
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.0520 
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10 


4 
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0. 
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Fig. 6. — Distribution of SN la occurrence as a 
function of peak /-band magnitude and redshift. 
Contours of predicted rates are drawn with heavy 
curves at 1 and 10, and light curves at intervening 
intervals of 2 and 5. The ridge line of an unred- 
dened "normal" SN la is drawn as a dashed line. 
The upper panel shows the distribution based on 
the SN la luminosity function if there were no 
reddening; the middle panel shows the effect of 
convolving with the model extinction distribution; 
and the bottom panel shows the expected distri- 
bution for our search parameters. 
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Fig. 7. — Cumulative number of SN la discovery 
is shown as a function of redshift. The histogram 
is the observed SN la; the curve is the prediction 
of what we should have seen given a flux limit of 
/ < 24 mag and a 31 d difference search. We think 
that our search is nearly complete for z < 0.6, but 
do not claim completeness for z > 0.6, because of 
the necessity to pick a subset for spectroscopy. 
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Table 15 — Continued 
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z 
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log( 


cz) 
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70 
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2 
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68 
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5 
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85 
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3 
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59 
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4 
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54 
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3 
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3.149 


0.030 





00 


R 


sn94ae 


225 
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59 
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25 


5 
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P 
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563 





3720 
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25 


5 


047 
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0.057 






P 
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69 


409 


-49 


082 





3780 
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25 


5 


054 


5.195 


0.081 






P 
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230 


027 


39 


659 





0077 


N2962 


30 


3 


363 


3.373 
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994 


30 
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N2441 


30 


3 
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2 


18 


RS 
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259 


953 


43 
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28 


5 
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02 


RT 
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28 
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4 
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58 
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30 
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76 
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-56 
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4 
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R 
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3 
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50 
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30 


3 


255 
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29 


4 
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29 
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5 
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5 
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22 
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5 
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99 
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65 
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5 
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62 
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5 
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04 
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60 
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28 


5 
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5 
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5 
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48 
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5 
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48 
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25 


5 
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5 
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11 
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59 
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3 
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40 
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4 
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16 
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56 
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3 
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3 


11 


s 
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69 


240 


36 


620 
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Anon 


28 


5 


120 


5.228 


0.041 





03 


RT 
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54 


608 
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28 


5 
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5.323 


0.047 





05 


RT 


sn97ck 


57 
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38 


452 
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28 


5 


464 


5.638 
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17 


RT 
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9 
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69 
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3 
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0.026 





06 


JT 
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8 


4 
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4.792 
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45 


S 
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265 
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8 


4 


269 


4.240 


0.045 
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S 
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-49 
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11 


3 
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3.674 
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04 
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-33 
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Anon 


11 


3 


950 


4.033 
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25 
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25 
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11 


3 
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32 
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8 


4 
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S 
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87 
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11 


3 
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0.043 





46 


JT 
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377 
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8 


4 
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4.254 
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71 


S 
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8 


4 
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0.039 





10 


S 


sn97fd 
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8 


4 
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0.065 





00 


s 


sn97ff 


125 
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54 
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1 


7550 
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32 


5 


721 


5.905 


0.083 





00 


R 


sn98D 


63 


777 


72 


905 





0132 


N5440 


11 


3 


597 


3.605 


0.069 





20 


J 


sn98I 


216 


646 


18 


718 





8870 


Anon 


34 


5 


425 


5.575 


0.053 





25 


JST 


sn98J 


238 


498 


32 


106 





8330 


Anon 


34 


5 


397 


5.615 


0.052 





26 


JST 


sn98M 


260 


707 


60 


423 





6300 


Anon 


34 


5 


276 


5.396 


0.045 





28 


RST 


sn98V 


43 


942 


13 


346 





0170 


N6627 


11 


3 


707 


3.686 


0.028 





31 


JST 


sn98ab 


124 


861 


75 


194 
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N4704 


11 


3 


921 


3.876 


0.025 





59 


JST 


sn98ac 


239 


012 


31 


832 





4670 


Anon 


34 


5 


146 


5.260 


0.043 





27 


RST 


sn98aj 


238 


163 


31 


312 





8600 
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34 


5 


411 


5.698 


0.070 





20 


ST 


sn98aq 


138 


837 


60 


269 





0037 


N3982 


1 


3 


045 


3.181 


0.023 





06 


RST 


sn98bp 


43 
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20 


481 





0104 


N6495 


11 


3 


494 


3.475 


0.029 





42 


RJST 


sn98br 


6 
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50 


757 
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A2029 


8 


4 


385 


4.442 


0.063 





39 


S 
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Table 15 — Continued 



SN 


l" 


b n 


z 


Host 


Ref. a 


log(cz) 


(log(dffo)) 




± 


(Ay) 


Contrib 


sn98bu 


234.413 


57. 


,020 


0.0043 


N3368 


12 
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2.912 





.028 
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RST 


sn98cm 


332.822 


62. 


.779 


0.0800 


A1780 


8 


4.380 


4.398 





.059 


0.36 


S 


sn98co 
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-44. 


.941 


0.0171 
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11 
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.046 


0.28 


JT 
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11 
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0.03 


RT 
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0.0157 
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3.727 





.024 


0.68 


JST 


sn98dh 


82.828 


-50. 


.644 


0.0076 


N7541 


11 


3.358 


3.403 





.034 


0.50 


ST 


sn98dk 


102.856 


-62. 


.161 


0.0120 


U00139 


11 


3.556 


3.571 





.028 


0.80 


JST 


sn98dm 


145.975 


-67. 


.406 


0.0054 


M+010444 


11 


3.209 


3.498 





.028 


1.06 


ST 


sn98do 


135.234 


-61 


,849 


0.0920 


Anon 


8 


4.441 


4.462 





,069 


0.00 


S 


sn98dv 


272.283 


-40. 


.300 


0.1550 


Anon 


8 


4.667 


4.706 





,124 


0.42 


S 


sn98dw 


143.355 
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0.0490 


A0151 


8 


4.167 


4.208 





.051 


0.00 


S 


sn98dx 


77.674 


26. 


,667 


0.0538 


U11149 


11 


4.208 


4.175 





.030 


0.05 


JT 


sn98dz 


247.207 


-56. 


.375 


0.0910 
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8 
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4.458 





.045 


0.16 


S 
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0.0570 


A3266 


8 
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4.314 





.043 


0.58 


S 
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.711 


0.0200 
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11 
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.039 


0.57 
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sn98ef 
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.566 


0.0170 


U00646 


11 
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.028 


0.08 
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sn98eg 
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-42. 
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0.0234 
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11 
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.042 


0.26 


JT 
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143.189 
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S 
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0.02 


ST 


sn99X 


186.585 
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11 
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3.888 





,037 
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JT 
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30 
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0.0157 


N2595 
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0.0099 
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A3392 


8 
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4.242 





.039 
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S 


sn99by 
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0.0027 


N2841 
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2.908 


3.067 





.167 


0.72 


ST 


sn99cc 


59.667 


48. 


.740 


0.0316 


N6038 


11 


3.977 


3.983 





.029 


0.13 


RJST 


sn99cl 


282.270 


76. 


.509 


0.0082 


N4501 


15 


3.391 


3.009 





.090 


1.68 


RS 


sn99cp 


334.851 


52. 


.708 


0.0104 


N5468 


15 


3.494 


3.525 





,026 


0.04 


RJT 


sn99cw 


101.768 


-67. 


,906 


0.0112 


M-010201 


11 


3.526 


3.469 





.034 


0.17 


J 


sn99da 


89.732 


32 


,649 


0.0121 


N6411 


14:20 


3.560 


3.623 





.044 


0.58 


RJST 


sn99dk 


137.349 


-47. 


.464 


0.0141 


U01087 


14 


3.626 


3.686 





.039 


0.07 


RT 


sn99dq 


152.839 


-35 


,871 


0.0136 


N0976 


11:20 


3.610 


3.548 





.023 


0.35 


JST 


sn99ef 


125.719 


-50. 


.086 


0.0380 


U00607 


11 


4.057 


4.147 





.033 


0.03 


JT 


sn99ej 


130.442 


-28. 


.946 


0.0127 


N0495 


11 


3.581 


3.684 





.031 


0.08 


JT 


sn99ek 


189.403 


-8. 


.234 


0.0176 


U03329 


11 


3.722 


3.687 





.044 


0.62 


JT 


sn99ff 


168.784 


-52. 


.949 


0.4550 


Anon 


35 


5.135 


5.289 





.048 


0.16 


RST 


sn99fh 


166.731 


-53. 


.793 


0.3690 


Anon 


35 


5.044 


5.135 





.061 


0.70 


ST 


sn99fj 


166.879 


-53. 


.734 


0.8150 


Anon 


35 


5.388 


5.564 





.051 


0.09 


RST 


sn99fk 


166.405 


-53 


,167 


1.0560 


Anon 


35 


5.500 


5.659 





,054 


0.03 


RST 


sn99fm 


167.075 


-52 


,994 


0.9490 


Anon 


35 


5.454 


5.580 





.053 


0.04 


RST 


sn99fn 


188.216 


-31 


,901 


0.4770 


Anon 


35 


5.155 


5.249 





.039 


0.17 


RST 
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4. Cosmological Implications 
4.1. Hubble Diagrams 

At a redshift of 0.5, the difference in luminos- 
ity distance between a flat and open universe with 
n M = 0.3 is 0.22 mag. With 230 SN la and 79 
at redshifts greater than 0.3, we might expect to 
be able to distinguish between these cosmologies 
at better than 6tr. The real concern, however, is 
whether we are making a systematic error or hit- 
ting a systematic noise floor for SN la, and that 
we are fooling ourselves in thinking that we are 
making progress by averaging together more and 
more SN la. 

The goal of this particular campaign was to 
strive to find and follow SN la at distinctly higher 
redshift than the majority of results at z = 0.5, 
since the effect of a dark-energy dominated cos- 
mology looks very different from most plausible 
systematic effects at z « 1. If the universe really 
has no dark energy and the dimness of SN la by 
0.22 mag at z = 0.5 is the result of a systematic 
error proportional to the redshift, one might ex- 
pect SN la to be dimmer by 0.44 mag at z = 1 
instead of 0.26 mag which would be the case for a 
flat universe. 

Figure 8 shows the results from this campaign 
plotted over the full set of data from Table 15. 
Evidently there is a tendency for the SN la to be 
brighter at z > 0.9 (although the uncertainties are 
still large), suggesting that we are truly seeing a 
change in the deceleration of the universe and the 
effects of dark energy. 

Given the heterogeneous nature of the data and 
the continual worry about non-Gaussian distribu- 
tions, it is difficult to be quantitative about how 
well we can reject the possibility that we are seeing 
a systematic error in SN la brightness in a universe 
with no dark energy, and we discuss this in more 
detail below. 

The scatter of the points from Table 15 obscures 
the trend with redshift, so we take medians over 
redshift bins and illustrate the result in Figure 9. 
Our redshift bins are derived from the data by 
choosing bins which have at least 12 SN la and a 
bin width of more than 0.25 in log z but no more 
than 50 SN la. The uncertainties, estimated from 
the 68% percentiles of the scatter of the points in 
each bin, are all less than 0.05 mag except for 0.08 




7, 



Fig. 8. — The Fall 1999 and other data points are 
shown in a residual Hubble diagram with respect 
to an empty universe. From top to bottom the 
curves show (^m^a) = (0.3,0.7), (0.3,0.0), and 
(1.0,0.0), respectively, and the Fall 1999 points are 
highlighted. 

mag in the highest z bin. We believe that this 
is approximately where systematics will start to 
dominate statistical error. That the points mostly 
lie below the (0.3, 0.7) cosmology gets at the heart 
of what we believe may the systematic uncertainty 
floor of the SN la distances in Table 15; we dis- 
cuss this below. Nevertheless, we are apparently 
starting to see SN la become brighter at z w 0.9, 
rather than becoming ever dimmer because of a 
systematic effect. 

4.2. Cosmological Parameters 

We can fit the data from Table 15 with a vari- 
ety of cosmological models. We do not want to fit 
all 230 SN la because the nearby ones are affected 
by peculiar motions which are beyond the scope of 
this paper to characterize. However, we do need to 
allow for the velocity uncertainty in modeling, so 
when we compute x 2 we & dd 500 km s -1 divided 
by the redshift in quadrature to the distance error. 
For all fits we report x 2 at the minimum as a func- 
tion of H , since this is equivalent to marginalizing 
over H n . 

Fitting all 230 SN la yields X 2 = 252. If we 
restrict our fits to z > 0.01, we obtain x 2 — 206 
for 195 points. We are also very wary of "extin- 
guished" SN la since we are not certain we un- 
derstand host extinction and reddening very well, 
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Table 15 — Continued 
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JT 
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23.736 
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3.356 


0.023 
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-19.504 


0.0260 
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3.892 


3.933 


0.028 


0.15 
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166.353 


22.791 


0.0193 


N2320 


11 


3.762 


3.740 


0.039 


0.38 


JST 


snOObk 


295.292 


55.233 


0.0266 


N4520 


14 


3.902 


3.911 


0.030 


0.26 


RJT 
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149.096 


32.004 


0.0165 
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3.705 


0.030 


1.02 
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99.883 


42.165 


0.0360 


M+111925 


11 


4.033 


4.112 


0.032 


0.13 


RJST 


snOOcn 


53.445 


23.318 


0.0233 


1X11064 


11 


3.844 


3.847 


0.023 


0.27 
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136.506 
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0.0070 


N0524 
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3.320 


0.023 


0.02 


RST 
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126.834 
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0.0164 
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11 
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3.677 


0.023 


0.01 


JST 


snOOdz 
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-56.390 


0.5000 


Anon 


13 
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5.352 


0.053 


0.07 


RST 
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0.4200 
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0.88 
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0.4700 
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13 


5.149 


5.340 


0.043 
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RST 
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13 
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0.041 
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5.236 
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188.297 
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0.4900 


Anon 
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Fig. 9. — The Fall 1999 and other data points are 
shown in a residual Hubble diagram with respect 
to an empty universe. In this plot the highlighted 
points correspond to median values in eight red- 
shift bins. From top to bottom the curves show 
(Sl M ,Sl A ) = (0.3,0.7), (0.3,0.0), and (1.0,0.0), re- 
spectively. 

19 and we are mindful of the covariance between 
Ay and d. If we further restrict the sample to 
A v < 0.5 mag we get X 2 = 169.5 for 172 SN la, 
and we regard this as the most useful sample for 
deriving cosmological parameters. There are fur- 
ther cuts which could be made: exclusion of the 
SCP data points (where there is no information 
on host galaxy extinction) gives us x 2 = 121.6 
for 130 SN la. However the constraints on cos- 
mological parameters do not change significantly 
and we are concerned about acquiring systematic 
problems along with spurious precision. 

'Phillips et al. (1999) give in their Table 2 the Galactic 
and host reddening E(B — V) for nearby supernovae. If, 
as suggested by the 44 SN la studied by Jha (2002, §3), 
most SN la obey Lira's (1995) relation, then the colors at 
maximum and in the tail of the B — V color curve can give 
us a good handle on the reddening at optical wavelengths. 
However, while most galaxies may have dust like that in our 
Galaxy, with R v = A V /E(B - V) » 3.1 (Cardelli et al. 
1989), there are heavily reddened objects such as SN 1999cl 
whose host galaxy has dust with Ry < 2 (Krisciunas et al. 
2001). Furthermore, high-rcdshift supernovae are typically 
too faint to detect with sufficient precision 30 to 90 rest- 
frame days after maximum light. Finally, while a combi- 
nation of optical and infrared data gives us a good handle 
on the extinction A v (Krisciunas et al. 2000, 2001, 2003), 
SN la with z > 0.10 are too faint in the rest-frame //-band 
and AT-band even for 10-m class telescopes to apply the V 
minus infrared method of determining extinction. 



In the analyses that follow we use various sub- 
samples of the 172 SN la with z > 0.01 and A v < 
0.5 mag. For each calculation we determine x 2 as 
a function of Ho, CIm, and Given that x 2 /N 
is close to unity, we convert this to a probability 
as exp(— 0.5x 2 ), normalize, marginalize over Ho, 
and locate contours of constant probability density 
which enclose 68%, 95%, and 99.5% of the proba- 
bility. In those calculations where we additionally 
constrain cosmology by using the galaxy cluster- 
ing shape measurements, we use observations from 
the 2dF redshift survey (Percival et al. 2001), 
Q M h = 0.20±0.03, and H = 72±8 km s" 1 Mpc" 1 
(Freedman et al. 2001). In all cases, contours are 
calculated for a cosmological constant equation of 
state (w = —1). 

Figure 10 shows these contours for a sample of 
94 SN la (x 2 = 103) which include only the recent 
results from Leibundgut et al. (2003), Suntzeff et 
al. (2003), Jha et al. (2003b), and the data pre- 
sented here. The error ellipses are virtually iden- 
tical to those drawn by Riess et al. (1998) and 
Perlmutter et al. (1999), although the data and 
analyses are completely different, confirming the 
conclusion that the universe is accelerating. This 
independent sample gives a consistent result: the 
acceleration inferred several years ago was not a 
statistical fluke. 

Figure 11 shows these contours for the sample 
of 130 SN la (x 2 = 122) which includes all the 
SN la which have been published by the HZT. 

Figure 12 shows these contours for the sample 
of 172 SN la (x 2 = 170) which includes all the 
SN la for which we have either data or published 
results. 

In all cases, the derived probability contours 
with a 2dF prior are remarkably tightly con- 
strained to the fltot = 1 line which is being delin- 
eated so well by CMB observations. 

As with previous results from SN la, the ap- 
proximate combination J7a — is relatively well 
constrained. For these data 

Q A - 1.40m = 0.35 ±0.14. 

However, the presence of these higher-redshift ob- 
servations from Fall 1999 as well as SN 1997ff 
(Riess et al. 2001) has started to close off the 
contours from running away to very large fi. If 
we assume a flat universe as well as w = —1, the 
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Fig. 10. — Probability contours for fl\ versus VI m 
are shown at la, 2c, and 3a with to = — 1. We 
also give la, 2cr, and 3a contours when we adopt 
a prior of Hm/i = 0.20 ± 0.03 from the 2dF survey 
(Percival et al. 2001). These constraints use only 
the 26 new SN la at z > 0.3 (which are completely 
independent of any which have been used before 
for cosmological constraints). 

SN la data require VIm = 0.28 ±0.05, independent 
of large-scale structure estimates. 

We can also fit for the equation of state of the 
dark energy, p — wp, where w = 1/3 for radia- 
tion, w = for pressureless matter (e.g., cold dark 
matter), w = — 1 for a cosmological constant, etc. 
Restricting ourselves to fi tot = 1, Figure 13 shows 
the contours for fl M and w, with and without the 
2dF prior. 

Dark energy candidates with w > —2/3 are 
strongly disfavored by these data. Figure 14 shows 
the data from Figure 13 marginalized over Qm- 
The combination of SN la and constraints on £1 m 
from large scale structure puts a sharp constraint 
on how positive w can be. The 95% confidence 
limits on w lie in the range —1.48 < w < —0.72. 
If we additionally adopt a prior that w > —1, wc 
find that w < —0.73 at 95% confidence. These 
constraints are again very similar to recent results 
reported using the WMAP satellite: w < —0.78 at 
95% confidence for priors of w > —1 and VLm from 
the 2dF redshift survey (Spergel et al. 2003). 



Fig. 11. — Probability contours for fi\ versus VI m 
are shown at la, 2a, and 3a with w = — 1. We 
also give la, 2a, and 3a contours when we adopt 
a prior of Hm/i = 0.20 ± 0.03 from the 2dF survey 
(Percival et al. 2001). These constraints use the 
113 SN la published by the HZT. 

Since the gradient of Hq to is nearly perpen- 
dicular to the narrow dimension of the £1a-£1m 
contours, we obtain a surprisingly good estimate 
of Ho to, illustrated in Figure 15. We find that 
H to = 0.96 ± 0.04. Although independent mea- 
surements of these parameters do not approach 
this precision yet, for comparison we note that 
h = 0.72 ±0.08 from Frccdman et al. (2001) corre- 
sponds to a Hubble time of H^ 1 = 13.6 ± 1.5 Gyr. 
Krauss & Chaboyer (2002) have recently found a 
median age of 12.5 Gyr for 17 metal-poor globular 
clusters. If these globular clusters formed 0.6 Gyr 
after the Big Bang, then the universe has to = 13.1 
Gyr, consistent with a world model with h — 0.72, 
VIm = 0.3, and Vl\ = 0.7. The assumed globular 
cluster incubation time of 0.6 Gyr corresponds to 
an epoch of star formation we would observe now 
at z ss 8. Adopting an uncertainty of ± 2 Gyr 
for t , the product H t = 0.96 ± 0.19. We also 
note that remarkable agreement with the WMAP 
result for the age of the universe: 13.7 ± 0.2 Gyr 
(Spergel et al. 2003). 
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Fig. 12. — Probability contours for fl\ versus VI m 
are shown at la, 2c, and 3a with to = — 1. We 
also give la, 2a, and 3a contours when we adopt 
a prior of ClMh = 0.20 ± 0.03 from the 2dF survey 
(Percival et al. 2001). These constraints use the 
full sample of 172 SN la with z > 0.01 and A v < 
0.5 mag. 

5. Discussion 

5.1. Testing systematics 

Drell, Loredo, & Wasserman (2000) (DLW) 
wrote a very thoughtful critique of the then ex- 
tant SN la data, questioning whether a cosmolog- 
ical constant was required or a systematic error 
could account for the observations. We present 
here more than twice as many objects with im- 
proved error estimates. We believe that some of 
the points noted by DLW are simply caused by se- 
lection effects, for example the observed diminu- 
tion of the range of A with redshift. However, 
their point that a large systematic shift in the lu- 
minosity of supernovae with redshift, coupled with 
a large positive CIm could fit the data nearly as 
well as a universe with a cosmological constant, is 
still valid." 

Using our preferred sample of 172 SN la we 
reevaluate the Bayesian odds factor between a 
model "FPW" with no systematic error and a cos- 
mological constant and a "Model II" with no cos- 
mological constant but a systematic error which 



Fig. 13. — Probability contours for dark energy 
parameter w versus tt m are shown at la, 2a, and 
3cr when fi tot = 1. We also give la, 2a, and 3a 
contours when we adopt a prior of Om/i = 0.20 ± 
0.03 from the 2dF survey (Percival et al. 2001). 
This sample includes all 172 SN la with z > 0.01 
and Ay < 0.5 mag. 

goes as a power law of (1 + z) (i.e., the distance 
modulus has a systematic error term of /31n(l + z) 
- this is the DLW "Model II"). When we inte- 
grate both models over H , £Im, and either £1^ or 
(3 and divide by plausible prior ranges for these pa- 
rameters, we get an odds factor of 29 in favor of the 
FPW model. Most would interpret this as "pretty 
good" odds that the FRW model is to be preferred 
over a systematic error, but not extremely strong 
odds. 

However, in this Bayesian framework, the ra- 
tio of the probabilities for the two models also in- 
cludes the ratio of the prior probabilities. The 
systematic model can match the observations rea- 
sonably well, but it does so only for VIm ~ 1 and a 
(3 which corresponds to a systematic which rises to 
0.7 mag at z = 1. Both of these parameters seem 
out of step with other lines of investigation and 
would be quite disfavored by reasonable priors. 

We certainly do not expect that Om could be 
anywhere near unity. The evidence against this 
comes from the observations that the H to prod- 
uct appears to be much larger than 2/3, and there 
is ample evidence on many scales that there is 
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Fig. 14. — The probability distribution for w using 
the 2dF prior, marginalized over Qm- The 95% 
confidence limits are —1.48 < w < —0.72. 

nowhere near as much mass as VLm = 1 in clus- 
tered form (e.g., Blakeslee et al. 1999; Percival et 
al. 2001). 

The argument against such a very large sys- 
tematic error (0.7 mag at z = 1) is theoretical and 
empirical. We work with SN la precisely because 
they are dead stars which should have very little 
difference between today and 7 Gyr ago, and we 
certainly do not have any theoretical expectation 
that explosions have become a factor of 2 dimmer 
in the mean over the previous 7 Gyr. 

Empirically, we argue that the SN la we sec lo- 
cally should have at least the range of properties 
seen at z = 1 because there are low metallicity sys- 
tems today and there appear to be SN la which 
are exploding quite promptly (e.g., NGC 5253 at 
0.1 L* had a starburst recently and has had 2 
SN la in the past century). Likewise, if SN la 
were as much as a factor of 2 dimmer in the past, 
we would expect there to be a significant disper- 
sion — surely some of the SN la at z = 1 would 
produce explosions comparable to those of today. 

We are always cognizant of the possibility of 
observational systematic error or problems with 
our analysis, but 0.7 mag is considerably beyond 
what we think might be possible. 

The exact factor by which Model II would 
therefore be disfavored according to these prior 
considerations is subjective since the place where 
the systematic probability is becoming sizable is 
out on the tails of the priors we assign for Am and 
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Fig. 15. — The probability distribution for Hq to 
given the SN la observations is tightly constrained 
to 0.96 ± 0.04, and an approximating Gaussian 
curve. 

(3. We think it is not worthwhile to try to assign 
a quantitative probability to these, other than to 
note that the factor by which the net probability 
of Model II relative to Model FRW will be disfa- 
vored will be significantly higher than the factor 
of 29 from the integrated likelihoods. 

Therefore we believe that we have very strong 
evidence that a systematic effect which goes as a 
power law of (1+z) is not likely to match the SN la 
data. In particular, a simple gray or intergalactic 
dust model does not appear to be viable, and evo- 
lutionary effects such as metallicity, C/O ratio, or 
host dust properties which change smoothly and 
monotonically with z are also likely not to fit the 
data as well as an FRW cosmology. 

It is always possible to imagine a systematic 
error which is tuned to match the FRW turnover 
near z = 1, but it might be less a priori plausible 
than a simple (1 + z) power. However, a proper 
Bayesian analysis would have to take into account 
the volume of the space from which such a model 
was drawn and the model might be disfavored for 
that reason, even though it might fit the SN la 
data very well. 

5.2. Selection and Other Effects 

The Fall 1999 survey went deeper and redder 
than any previous campaign. As a result we did 
start to see extinguished SN la at z w 0.5 and we 
expect that we are less subject to selection bias. 
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Our model for rates seems to work remarkably 
well, suggesting that there is not a rare population 
of luminous objects which we do not know about 
from local samples. The exceptions are SN 1999fo 
(Gertie) and SN 1999fu (Alvin), which were not 
SN la, but were definitely unusual, flux-variable 
objects. 

At face value the average absorption of the 
distant supernovae is significantly smaller than 
for the nearby sample. Out to z < 0.3 the 
f 51 supernovae suffered an average absorption of 
(Ay)=0.35±0.50 mag, while the 37 distant super- 
novae with z > 0.3 average (j4y)=0.20±0.23 mag. 
There is a large scatter in these averages, but the 
offset is 0.15 mag. The distribution of absorp- 
tion values for the nearby supernova sample has 
a long tail to higher absorptions, which is miss- 
ing for the distant sample. However, the distri- 
bution of absorptions of the distant supernovae 
is much narrower. Most probably this is a sig- 
nature of the magnitude-limited searches, where 
heavily obscured supernovae are not detected. Se- 
lecting the 130 objects (172 less 42 SCP SN la 
without host extinctions) which were used to con- 
struct the likelihood functions (z > 0.01 and 
Ay < 0.5 mag), we find that the averages decrease 
to (Ay) = 0.17 ± 0.14 mag for the nearby super- 
novae and (Ay) = 0.14± 0.11 mag for the distant 
ones. The difference has decreased to less than 
0.03 mag and clearly the faintness of the distant 
supernovae cannot be attributed to errors in the 
treatment of absorption. In fact, more absorption 
would make the distant supernovae more luminous 
and move them to larger distances. The small av- 
erage absorption argues that the distant objects 
indeed are at larger distances than for an empty 
universe model, in agreement with the results from 
Riess et al. (2000), Leibundgut (2001), and Sulli- 
van et al. (2003). 

Of the objects at z > 0.8 in our large data set, 
several objects (principally those discovered be- 
fore 1999) rely on rest-frame {/-band data to mea- 
sure extinction. The rest-frame U -band data have 
a number of uncertainties, including less certain 
K-corrections, greater susceptibility to extinction, 
and less well studied behavior. Most importantly, 
the SN la have less intrinsic uniformity in the U 
bandpass. For the data set presented here, the 
average measured distance to those objects using 
the U band as a primary source of information for 



distance measurement is further than to those ob- 
jects which do not rely on [/-band data, but the 
difference is not statistically significant. These ef- 
fects will be further discussed in Jha et al. (2003c) 
and Suntzeff et al. (2003). 

We have not corrected these data for gravita- 
tional lensing. It should be negligible at these red- 
shifts, but will become much more serious as we 
push to higher redshift. 

Nevertheless, we do think we can start to per- 
ceive systematic errors in the data at the 0.04 
mag level. The contours seen in Figures 10-12 
have their most probable £Im, well offset from 
the f^M indicated by the 2dF survey and a flat 
universe. This is not a failure of our probabil- 
ity analysis. We have performed extensive Monte 
Carlo tests on simulated data sets and we recover 
Q.M, with very little bias. What drives our best 
fit solution to high values of Qm , is the fact 
that the data make a very fast transition from re- 
cent acceleration to earlier deceleration, i.e. the 
transition from the faintness at z ~ 0.5 to the 
brightness at z ~ 1 is very fast. The contours in 
Figures 10-12 are quite shallow in the Qm + ^a 
direction, and a systematic error of 0.04 mag is 
quite sufficient to move our most probable values 
for Qm,&a to these high values. This may be a 
statistical fluke at the 1.5-er level, although the 
systematic offsets of the points in Figure 9 from a 
cosmological line may point to systematic errors. 

We have described some of the possible culprits 
for such a systematic above. Selection biases be- 
cause of UV selection or any of the other practical 
difficulties in earring out a supernova search are 
always a worry, particularly since they are so hard 
to quantify and because SN la are so non-uniform 
in the UV. Photometric and K-correction errors 
are certainly possible, since absolute photometry 
of objects with non-stellar SEDs at m > 20 is ex- 
tremely hard to achieve at the 0.01 mag level. It is 
conceivable that gravitational lensing is more im- 
portant than has been estimated; we expect that 
the SN la found by the GOODS survey will offer 
a strong test of this. Most likely there is a vari- 
ety of errors throughout this very heterogeneous 
set of data, and the best way forward is acquisi- 
tion of new data from well designed observations 
which are taken and analyzed with more care. Al- 
though the basic conclusion that the Universe is 
accelerating is not in jeopardy since the Om — ^a 
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direction is relatively insensitive to systematic er- 
ror, we remain wary of over interpretation of the 
extant SN la data and urge others to do the same. 

5.3. Comparison with the CMB 

As this paper was nearing completion, the 
WMAP team reported new results from a high an- 
gular resolution all-sky map of structure (Spergel 
et al. 2003). Their main conclusions, namely 
that we live in a geometrically flat Universe with 
Mtotai, Am, and fi A of 1.00, 0.73, and 0.27, re- 
spectively, are entirely consistent with the results 
presented here. This convergence of results lends 
great strength to what is emerging as the new 
"Standard Model" of contemporary cosmology, 
in which the role of matter is secondary to the 
accelerating expansion driven by dark energy. 

Ascertaining the physical mechanism that is re- 
sponsible for this accelerating expansion is one of 
the crucial next steps for observational cosmology. 
Since supernovae are well suited to probing pre- 
cisely the redshift range over which A begins to 
dominate, new projects are being undertaken to 
exploit this. 

5.4. What Next? The Future 

Future campaigns to observe SN la at high red- 
shift will focus on three themes: refining measures 
of luminosity distance, understanding and control- 
ling systematic errors, and learning about SN la 
physics. The payoff for measuring luminosity dis- 
tances accurately is that we can determine the 
equation of state of the dark energy and whether it 
has changed as a function of time (redshift). There 
is certainly a systematic floor to our ability to de- 
termine distances from SN la, but we do not yet 
know what it is. It appears that gray dust is not 
an issue, but gravitational lensing, for example, is 
certainly a contaminant. 

Opinions differ as to where the incompletely 
understood SN la physics causes a cosmic scat- 
ter which does not average out, but most would 
agree it lies in the range of 0.02-0.04 mag. This 
implies that the luminosity distance at a given 
redshift is not improved by observing more than 
25-100 SN la at that redshift. There are many 
outstanding questions about SN la physics which 
may best be approached by observations of SN la 
at high redshift. Measuring SN la rates and com- 



parison with the star-formation rate may elucidate 
the SN la explosion mechanism, SN la progenitors, 
and gestation period. The distribution of SN la 
parameters as a function of redshift may also be 
useful for understanding SN la and host galaxy 
physics. 

The HZT, in collaboration with the Surf's-Up 
group at the IfA, undertook a six-month survey 
of 2.5 deg 2 in i?, I, and Z. Many SN la as well 
as other transient objects were discovered and 
followed (Barris ct al. 2003). Members of the 
HZT have begun the ESSENCE survey (Smith 
ct al. 2002) to discover and follow 200 SN la in 
the range 0.1 < z < 0.8 with the goal of mea- 
suring the equation of state parameter w of the 
dark energy to ±10%. Riess ct al. (2003) have 
an ongoing program with HST and the GOODS 
survey, using the Advanced Camera for Surveys 
(ACS) to discover and follow approximately 6 
SN la at z ~ 1.4. This will be exceedingly pow- 
erful for measuring the effects of dark energy, 
understanding systematics, and constraining any 
time evolution of w. The CFHT Legacy Survey 
(http : / / www. cfht . hawaii . edu /Science / CFHLS / ) 
has the goal of finding and following ~ 2000 
SN la in multiple colors. The Lick Observa- 
tory Supernova Search (LOSS; Filippenko et al. 
2001) is finding and monitoring record numbers 
of very nearby supernovae, with 82 in 2002 alone 
(http :/ /astro . berkeley. edu /<~bait /kait . ht ml) . The 
Nearby Supernova Factory (Aldering et al. 2002) 
will obtain spectrophotometry for ~ 300 nearby 
SN la, suitable for K-corrections of high-redshift 
supernovae and for exploring SN la physics and 
systematics. In the near future, Pan-STARRS 
(Kaiser et al. 2002) should discover and follow a 
vast number of SN la (>10 4 per year). Finally, in 
the more distant future, the proposed SNAP satel- 
lite (Nugent 2001) would follow several thousand 
objects with an experiment designed to minimize 
observational systematic errors, limited only by 
the floor imposed by nature. The future indeed 
looks bright for cosmological and related studies 
of SN la. 
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6. Appendix A 

6.1. N(N — l)/2 differences 

We use the Alard & Lupton (1998) code to 
convolve one image to match the other. We fit 
a suitable star in one of these images with the 
Schechter "Waussian" (inverse of a truncated Tay- 
lor series for e x ) to derive a standard PSF. Since 
we know where the supernova is located, we do a 
two-parameter fit (peak and sky) of this standard 
PSF, and we scale the flux measured in the stan- 
dard star by the ratio of the fit peaks to derive 
a supernova flux difference for these two images. 
We also select roughly ten other locations in the 
image and insert a copy of the standard star scaled 
to match the supernova flux, and pick it up again 
using the same procedure of two-parameter fits. 
This tells us of any bias (we do not see any) and 
gives us an estimate of the flux error for this pair 
of images. 

Applying this procedure to all N(N— 1)/2 pairs 
creates an N x N antisymmetric matrix of flux 
differences, and a symmetric matrix of errors. We 
also choose approximately a dozen stars for each 
supernova field which are bright enough to have 
good S/N in all images but faint enough not to be 
saturated, and perform a PSF flux measurement 
for all stars in all images. Matching up these tables 
of stars, we can form a mean flux scale between all 
pairs of images, and then assemble another N x 
N antisymmetric matrix of magnitude differences. 
The flux differences are all put on a common scale 
using the magnitude differences. 

Finally, we fit these antisymmetric matrices as 
the difference between the i th and j th entries of 
an N vector. We derive two sets of error esti- 
mates: the first is based on the pairwise flux error 
estimates and the second comes from the resid- 
uals between the observed difference matrix and 
the matrix assembled by differencing the terms of 



the resultant vector. These estimates are gener- 
ally consistent with one another, although the ob- 
served errors reflect the uncertainties of the indi- 
vidual observations and the mismatch errors give 
a broader view of the inconsistency of any obser- 
vation with the others. 
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